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ABSTRACT 


Exterior  scattering  from  a  three-dimensional  impedance  body 
can  be  formulated  in  terms  of  various  integral  equations  derived 
from  the  Leontovich  impedance  boundary  condition  (IBC).  The 
electric  and  magnetic  field  Integral  equations  are  ill-posed 
because  they  theoretically  admit  spurious  solutions  at  the 
frequencies  of  interior  perfect  conductor  cavity  resonances. 

A  combined  field  formulation  is  well-posed  because  it  does  not 
allow  the  spurious  solutions.  This  report  outlines  the 
derivation  of  IBC  integral  equations  and  describes  a  procedure 
for  constructing  moment-method  solutions  for  bodies  of 
revolution.  Numerical  results  for  scattering  from  Impedance 
spheres  are  presented  which  contrast  the  stability  and  accuracy 
of  solutions  to  the  ill-posed  equations  with  those  of  the 
well-posed  equation.  The  results  show  that  numerical  solutions 
for  exterior  scattering  to  the  electric  and  magnetic  field 
Integral  equations  can  be  severely  contaminated  by  spurious 
resonant  solutions  regardless  of  whether  the  surface  impedance  of 
the  body  is  lossy  or  lossless.  Solutions  to  the  combined  field 
equation  are  not  contaminated  by  the  spurious  solutions  when  a 
suitable  choice  of  Integral  equation  weighting  coefficient  can  be 


determined.  However,  the  dete 


S  nation  of  th 
ance  case  is 


the  weighting 


coefficient  for  the  general  Impedance  case  is  more  difficult  and 
more  critical  than  for  the  perfectly  conducting  scatterer  case. 
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INTRODUCTION 


A  well-known  method  for  modeling  electromagnetic 
Interactions  with  certain  types  of  Imperfectly  conducting  bodies 
Involves  describing  the  body  surface  with  a  Leontovlch  impedance 
boundary  condition  (IBC).1"4  This  approximation  makes  integral 
equation  formulations  of  scattering  problems  nearly  as  simple  as 
those  for  perfectly  conducting  bodies.5  Until  the  present, 
solutions  to  IBC  integral  equations  for  exterior  scattering  from 
generally-shaped,  three-dimensional  bodies  have  been  obtained 
from  either  the  magnetic  field  integral  equation  (MFIE)6"8  or  the 
electric  field  integral  equation  (EFIE).9  These  two  IBC 
equations  can  be  described  as  being  mathematically  ill-posed  for 
exterior  scattering  problems  because  at  certain  discrete 
frequencies,  the  integral  operators  of  each  become  singular  and 
the  equations  no  longer  have  unique  solutions.10  These  discrete 
frequencies  are  the  interior  resonances  of  a  perfectly  conducting 
cavity  having  the  same  shape  as  the  impedance  target.10 
Numerical  solutions  to  either  the  MFIE  or  EFIE  near  a  cavity 
resonance  will  result  in  ill-conditioned  matrices  and  scattering 
results  which  may  be  contaminated  by  spurious  solutions  due  to 
the  cavity  resonances.*  In  contrast  to  the  MFIE  and  EFIE,  the 

*  The  terminology  "spurious  solutions"  refers  to  unwanted 
mathematical  solutions  which  have  no  physical  connection  with  the 
exterior  scattering  problem. 
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IBC  combined  field  integral  equation  (CFIE)  is  mathematically 
well-posed  because  it  does  not  allow  the  spurious  resonant 
solutions. 1 0 

This  report  investigates  numerical  solutions  for  scattering 
from  an  impedance  body  of  revolution  (BOR)  as  determined  by  the 
MFIE,  EFIE  and  the  CFIE.  The  principal  objectives  are  to:  (1) 
describe  the  method  by  which  moment-method  solutions  to  the  three 
equations  can  be  constructed,  (2)  contrast  the  numerical 
stability  and  accuracy  of  solutions  to  the  MFIE  and  the  EFIE 
with  those  of  the  CFIE,  and  (3)  investigate  the  sensitivity  of 
the  CFIE  to  parameters  such  as  integral  equation  weighting 
coefficient  and  moment-method  segmentation  for  various  types  of 
surface  impedance  (i.e.  reactive,  lossy,  etc.).  The  report 
begins  with  a  brief  discussion  in  Section  2  of  the  derivation  of 
the  three  equations  and  the  well-posedness  of  each.*  In 
Section  3  it  will  be  described  how  existing  moment- 
method  numerical  algorithms  applicable  to  scattering  from 
perfectly  conducting  BORs11”13  can  easily  be  generalized  to  the 
IBC  case.  A  Fortran  IV  computer  program  which  contains  the 
modifications  discussed  in  Section  3  is  listed  in  Appendix  A. 
Numerical  results  which  address  objectives  (2)  and  (3)  above  are 


*  A  general,  theoretical  treatment  of  the  well-posedness  of 
IBC  Integral  equations  is  given  in  Ref.  10.  The  numerical 
results  of  this  report  represent  the  practical  implementation  of 
one  of  the  methods  discussed  in  Ref.  10. 
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presented  in  Section  4.  These  results  are  for  scattering  from 
impedance  spheres,  and  in  all  cases  the  moment-method  solutions 
were  compared  to  an  "exact"  series  solution.14 

Some  conclusions  concerning  the  results  are  presented  in 
Section  5.  Briefly,  the  results  demonstrate  that  numerical 
solutions  to  the  MFIE  and  EFIE  for  exterior  scattering  can  be 
severely  contaminated  by  spurious  solutions  regardless  of  whether 
the  surface  impedance  is  lossy  or  lossless.  The  errors  in  pre¬ 
dicting  scattering  from  bodies  having  non-zero  surface  impedances 
are,  in  general,  larger  than  for  the  perfectly  conducting  case.* 
Also,  the  results  demonstrate  that  the  solution  accuracy  of  the 
CFIE  as  a  function  of  integral  equation  weighting  coefficient  and 
moment-method  segmentation  is  strongly  influenced  by  the  value  of 
surface  impedance.  Criteria  established  for  the  perfect 
perfect  conductor  case7*11  may  not  work  well  when  the  surface 
impedance  is  non-zero. 


*  A  perfectly  conducting  surface  can  be  considered  as  the 
special  case  of  an  impedance  surface  when  the  surface  impedance 
equals  zero. 


2.0  REVIEW  OF  THE  FORMULATION  AND  STABILITY  OF  IMPEDANCE 

BOUNDARY  CONDITION  INTEGRAL  EQUATIONS 

The  physical  scattering  problem  is  depicted  in  Fig.  2.1.  A 
uniform  plane  electrc magnetic  wave  (fi^.fi*-)*  is  incident  upon  an 
imperfectly  conducting  body  surrounded  by  a  closed  surface  S. 
ce1.*1)  are  assumed  known  and  have  a  monochromatic  time  variation 
of  eJ“c.  When  the  incident  wave  interacts  with  the  body, 
volume  electric  current  density,  JVol»  induced  inside  S  and 
radiates  a  scattered  field  (£8,ft8).  Assuming  that  the  radiation 
conditions  are  met  and  the  target  is  passive,  the  scattered 
fields  are  unique. 

In  solving  the  problem  depicted  in  Fig.  2.1,  it  is  often 
useful  to  apply  the  Equivalence  Principle15  which  states  that  the 
body  with  volume  electric  currents  can  be  conceptually  replaced 
by  a  sourceless  volume  having  electric  (3)  and  magnetic  (ft) 
surface  current  densities  on  S.  The  jump  conditions  which  must 
be  satisfied  on  S  are:15 

n  x  (£+  -  £_)  -  -  ft  (2.1) 

n  x  (fi+  -  ft_)  -  J  (2.2) 

*  The  field  point  dependence  (r)  of  all  vectors  and  dyadics  is 
suppressed  in  most  cases  throughout  the  text. 
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where  (£t,fi±)  are  the  total  electric  and  magnetic  fields  just 

A 

outside  (+)  or  just  inside  (-)  the  surface  and  n  is  the  outward 
unit  normal  to  the  surface.  These  total  fields  are  the  sums  of 
incident  fields  due  to  sources  at  infinity  and  scattered  fields 

due  to  (3, ft): 


fi±  -  fi1  +  s® 

(2.3) 

a±  -  a1  +  a® 

(2.4) 

In  constructing  appropriate  integral  equations,  a  certain 
arbitrariness  exists  in  the  choices  of  (3, ft)  and  (fi_,fi_). 

Although  the  exterior  fields  must  be  unique,  there  are  infinitely 
many  sets  of  equivalent  currents  and  interior  fields  which  will 
give  rise  to  the  correct  exterior  fields.  Thus,  one  can  choose 

(£_,$_)  when  formulating  the  problem  in  order  to  obtain  unique 
equivalent  currents.  As  depicted  in  Fig.  2.2,  a  convenient 
choice  for  the  interior  fields  is  often  the  null  field. 

Assume  that  an  imperfectly  conducting  body  has  the  following 
Leontovich  impedance  boundary  condition  (IBC)  satisfied  on  its 
surface1 ~4 : 

n  x  x  n  ■  n  2g  •  (n  x  ft+)  (2.5) 

where  lg  is  the  normalized  dyadic  surface  impedance  of  the  body 
containing  information  about  the  material  properties  within  S  and 


n  is  the  free  space  wave  impedance.  Upon  substitution  of 
Eqs.  2.1  and  2.2  into  Eq.  2.5.  the  following  condition  relating 
the  tangential  components  of  the  interior  fields  to  the 
equivalent  currents  is  obtained: 

£_xn+$- 

-n  n  x  2s  •  ^  -  n  n  x  5s  .  (n  x  ft_)  (2.£) 

When  choosing  the  null  field  problem  depicted  in  Fig.  2.2,  the 
following  two  conditions  must  hold  on  the  surface: 

n  x  *  0  (2.7) 

n  x  “  0  (2.8) 

For  this  choice  of  (£_,$_)  a  simple  relationship  for  ft  in  terms 
of  3  is  obtained: 

ft  -  -n  n  x  Zs  •  J  (2.9) 

This  last  relation  effectively  reduces  the  problem  from  two 
coupled  vector  integral  equations  in  two  vector  unknowns,  (5, ft), 
to  one  vector  integral  equation  in  one  vector  unknown,  5.  Note 
that  each  vector  integral  equation  consists  of  two  coupled  scalar 
Integral  equations  in  two  scalar  unknowns.  Equations  2.7  and  2.8 


w 


each  lead  to  a  different  vector  integral  equation  for  the 
unique  source  3.  Using  equations  2.3  and  2.4  and  integral 
representations  for  and  fi®  in  terms  of  (J.ft),15*16 
Eqs.  2.7  and  2.8  lead  to  the  following  IBC  electric  field 
integral  equation  (EF1E)  and  IBC  magnetic  field  integral 
equation  (MFIE),  respectively: 


nxfci«-Y  +  nx/a  U“U  J(r')  ♦ 


+  ft(f')  x  V*  +  jlj  7’*]  ds' 


EFIE 

(2.10) 


nxft1-  7"nx/s  ft(r')  • 


+  J(r')  *  v'*  -  jt-  7'.  ft(*')  7 ' • ]  ds’ 


where 


and 


fi  •  -ii  n  x  !  <  3 


1  e~^k  l*~r  I 
•{  -f 


MFIE 

(2.11) 


Equations  2.10  and  2.11  are  valid  in  the  limit  as  r+S  from  the 
inside.  The  terminology  "EFIE"  and  "MFIE"  in  this  report  will 
denote  the  generalized  impedance  electric  and  magnetic  field 


9 


integral  equations.  The  perfect  conductor  equations  are  the 
special  cases  when  2g  =  0. 

By  the  choice  of  vanishing  interior  fields,  two  different 
vector  integral  equations  were  obtained  for  determining  a  unique 
source  J.  Although  it  would  appear  that  either  the  EFIE  or  MFIE 
could  be  solved  for  3,  the  solution  must  ensure  that  both 
equations  are  satisfied.  This  is  required  because  Eq.  2.9  (used 
in  the  solution  for  J)  is  valid  only  if  both  Eqs.  2.7  and  2.8 
are  satisfied.  Here-to-fore,  it  has  been  assumed  that  a  solution 
to  either  the  MFIE  or  the  EFIE  will  satisfy  both  equations. 
However,  that  assumption  is  valid  only  away  from  discrete 
frequencies  located  at  the  perfect  conductor  (p.c.)  interior 
cavity  resonant  frequencies.10  At  a  p.c.  cavity  resonant 
frequency,  the  right  hand  sides  of  both  Eqs.  2.10  and  2.11  are 
singular  allowing  non-unique  mathematical  solutions  to  exist  at 
that  frequency.10  Consequently,  sets  of  spurious  equivalent 
currents  which  are  solutions  to  the  homogeneous  perfect  conductor 
equations  can  contaminate  the  particular  scattering  solution. 
Thus,  Eqs.  2.10  and  2.11  are  not  equivalent  to  each  other  at  a 
p.c.  cavity  resonant  frequency  and  the  solution  technique  must 
ensure  that  both  are  satisfied. 

Although  uniqueness  problems  with  the  EFIE  and  the  MFIE 
occur  at  discrete  frequencies,  it  is  still  advantageous  to  have 
an  integral  equation  formulation  which  ensures  a  unique  solution 
at  all  frequencies.  This  is  due  to  the  fact  that  numerical 
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solutions  to  either  ill-posed  equation  will  be  contaminated 
within  some  bandwidth  about  a  spurious  resonance  (due  to  the 
approximation  of  an  integral  equation  with  a  matrix  equation, 
numerical  roundoff  error,  etc.).  Since  the  perfect  conductor 
cavity  resonances  of  a  generally-shaped  target  are  not  known  a 
priori  and  become  increasingly  closely  spaced  with  increasing 
frequency,  it  is  apparent  that  for  large  scattering  targets  (with 
respect  to  a  wavelength)  the  EFIE  and  the  MFIE  may  be  unreliable 
at  all  frequencies. 

If  Eqs.  2.7  and  2.8  lead  to  two  equations  for  a  unique 
source  J,  then  a  linear  combination  will  also  yield  a  solution: 

n  x  ft  +  f  •  n  x  ft  “0  (2.12) 

“  a  " 

where  2  is  some  arbitrary,  generalized  weighting  coefficient 
which  could  be  a  function  of  r.  When  substituting  Eqs.  2.3  and 
2.4  into  2.12,  along  with  the  appropriate  integral 
representations  for  ft8  and  ft8,  one  obtains  an  IBC  combined  field 
integral  equation  (CFIE).  It  can  be  shown  that  a  solution  to  the 
CFIE  will  ensure  that  both  2.7  and  2.8  are  satisfied  at  all 
frequencies  when  2a  is  not  reactive.10  Thus,  the  CFIE  for  the 
impedance  body  case  rejects  spurious  resonances  in  a  similar 
fashion  as  does  the  perfect  conductor  CFIE.11*12 
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3.0  GENERALIZATION  OF  PERFECT  CONDUCTOR  MOMENT-METHOD  ALGORITHMS 

TO  THE  IMPEDANCE  BODY  CASE 

3.1  IBC  Integral  Equations  In  Perfect  Conductor  Integral 
Operator  Form 

In  Section  2,  three  Impedance  boundary  condition  Integral 
equations  were  derived  and  the  well-posedness  of  each  was  briefly 
discussed.  The  EF1E  and  MFIE  formulations  are  Ill-posed  due  to 
non-uniqueness  at  Interior  resonant  frequencies  of  the  perfectly 
conducting  cavity.  The  CFIE,  on  the  other  hand,  is  well-posed  at 
all  frequencies  (provided  2a  is  non-reactlve) .  In  this 
section,  it  will  be  shown  how  numerical  solutions  to  each  of 
these  equations  can  be  constructed  from  moment-method  matrix 
operators  applicable  to  scattering  from  a  perfectly  conducting 
body.  A  main  computer  program  which  accomplishes  this  procedure 
for  bodies  of  revolution  is  listed  in  Appendix  A. 

The  particular  perfect  conductor  algorithms  to  be 
generalized  are  those  of  Mautz  and  Harrington.12*13  These 
moment-method  programs  were  chosen  because  of  their  wide 
availability  and  good  documentation,  and  because  they  contain  all 
the  required  matrix  operators  for  solving  the  IBC  equations. 
Although  the  programs  are  restricted  to  bodies  of  revolution,  the 
basic  procedure  in  extending  these  programs  to  the  impedance  body 
case  would  also  be  applicable  to  moment-method  programs  for  more 
generally-shaped  bodies.  It  should  be  noted  that  the  procedure 
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to  be  described  in  this  section  has  already  been  accomplished  for 
the  EFIE  by  Iskander,  et.  al.9* 

It  is  appropriate  to  first  recast  the  IBC  integral  equations 
into  operator  notation  applicable  to  Mautz  and  Harrington's 
equations.  Their  integral  equations  for  scattering  from  a 
perfectly  conducting  body  can  be  written  as:11*12 

~  n  x  £*  x  n  -  L-.C5)  (3.1) 

n  a 

n  x  ft1  «  1^(3)  (3.2) 

where  the  integral  representations  of  the  operators  Lg  and  Lg 
are: 

Lgtf)  -  n  X  n  X  Jg  [jay  J(r')  «  +  v'»3(r')  7'*]ds'  (3.3) 

Ljjtf)  -  |  -  n  x  /s  J(r')  x  v'*  ds '  (3. A) 

Equations  3.1  and  3.2  are  the  perfect  conductor  special  cases  of 
the  IBC  EFIE  and  MFIE  equations  with  the  EFIE  put  into 
similar  vector  direction  and  units  as  that  of  the  MFIE. 

In  terms  of  the  perfect  conductor  operators,  Lg  and  Lg, 
the  IBC  integral  equations  are: 

*  The  formulas  in  Ref.  9  should  be  used  with  caution.  In 
addition  to  the  dyadic  impedance  boundary  condition  definition 
being  unclear,  there  appears  to  be  a  minus  sign  error  in  that 
definition  which  is  carried  throughout  all  subsequent  equations. 


13 


1  n  x  fil  x  n  -  Lgtf)  +  n  x  L^)  EFIE  (3.5) 

n  x  ft*  -  Ljj(3)  +  n  x  LgC^)  MFIE  (3.6) 

nxft*  +  a-^nxft*xn«  LH(J)  +  a  Lg(3) 

+  n  x  L-(-)  +  a  n  x  L„A  CFIE  (3.7) 

E>  ii  n  n 

where  a  is  a  constant  coefficient  in  the  CFIE. 

3.2  Review  of  the  Moment-Method  Solution  for  Scattering 
from  a  Perfectly  Conducting  Body 

Mautz  and  Harrington11-13  solve  Eqs.  3.1  and  3.2  for  a 

body  of  revolution  (BOR)  using  the  geometry  in  Fig.  3.1.  The 

AAA 

coordinates  (n,#,t)  represent  a  right-handed,  orthogonal  system 
for  describing  vectors  on  the  BOR  surface.  The  unknown  electric 

A  A 

surface  current  density,  consisting  of  t  and  t  components,  is 
expanded  in  terms  of  a  Fourier  series  in  +  and  in  terms  of 
overlapping  triangular  expansion  functions,  f,  in  t: 

3(?)  -  JC(t,t)  t  +  J*(t,*)  ; 
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Fig.  3.1.  Right-handed,  orthogonal,  curvilinear  coordinate  system  (ft,  $,  t) 
on  the  surface  of  a  body  of  revolution. 
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where. 


..  N  K  .  . 

JC(t.*)  -  l  l  fk(t)  ejn* 

n-1  k-1  nK  K 

-  I  I  I*  f,(t)  ejn* 

n-1  k-1  K 


(3.8) 


To  determine  the  number  and  location  of  all  fk(t),  the  top  half 
of  the  BOR  profile  is  segmented  into  short  straight  line 
segments.  The  number  of  expansion  functions,  K,  is  approximately 
1/2  the  number  of  straight  line  segments.  For  a  visual 
description  of  the  triangle  functions  see  Ref.  12. 

Mautz  and  Harrington  define  vector  testing  functions,  ft,  as: 


ft( r)  -  Wt"(t , ♦)  t  +  W*(t,*)  + 

where 

i-  *  N  K  .  A 

Wt(t,*)  -  W*(t,*)  -  l  l  f,(t)  e’jm*  (3.9) 

m-1  i-1  1 


Upon  substitution  of  the  expansion  for  3  into  Eqs.  3.1  and  3.2 
and  taking  the  appropriate  inner  products  of  the  equations  with 
the  testing  functions,  one  obtains  a  series  of  independent  matrix 
equations  corresponding  to  the  various  Fourier  modes.  (The 
Fourier  series  inner  product  is  non-zero  only  for  m-n  due  to  the 
azimuthal  symmetry  of  the  problem.)  For  the  n^  Fourier  mode: 
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<V  -  tznl  [In] 


lu  ]  -  [Y  ]  [I  1 
1  nJ  1  nJ  1  nJ 


(3.10) 

(3.11) 


or,  expressing  the  vector  nature  of  the  equations: 


tt 

ZC* 

n 

n 

♦t 

z** 

n 

n 

tt 

n 

n 

n 

n 

3.11 

are 

(3.12) 


(3.13) 


of  Eqs.  3.1  and  3.2  respectively.  The  known  vectors,  [Vn]  and 
[Un] ,  each  contain  2K  elements  which  are  obtained  from  the 
inner  products  of  the  incident  fields  with  the  testing 
functions.  The  vector  [In]  contains  the  2K  unknown  expansion 
coefficients  for  the  nc^  Fourier  mode  of  5.  [Zn]  and  [Yn] 
are  2K  x  2K  square  matrices  and  can  be  interpreted  as  normalized 
impedance  and  admittance  matrices  (not  to  be  confused  with 
surface  impedance  or  admittance).  Explicit  expressions  for  the 
elements  of  [Vn],  [UnJ,  [Zn]  and  [Yn]  are  given  in  Ref.  12.  For 
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each  Fourier  mode  retained  in  the  solution,  equation  3.10  or  3.11 
(or  Mautz  and  Harrington's  perfect  conductor  combined  field 
equation)  can  be  solved  for  the  [Inl  by  using  matrix  inversion 
or  factorization.  The  scattered  field  due  to  the  equivalent 
currents  of  each  mode  can  then  be  calculated  and  summed  up  for 
the  total  scattering  solution. 

3.3  Generalizing  the  Moment-Method  Solution  to  the 
Impedance  Body  Case 

For  the  IBC  case,  the  appropriate  integral  equations  are 
given  in  operator  form  in  Eqs.  3.5,  3.6  and  3.7.  These  integral 
equations  were  derived  from  the  impedance  boundary  condition  of 
Eq.  2.5  assuming  vanishing  interior  fields.  The  derivations  led 

to  an  expression  for  ft  in  terras  of  5,  restated  here  for 
convenience: 

ft  -  -n  n  x  Zs*  3  (2.9) 

Assume  that  the  normalized  dyadic  surface  impedance  is  diagonal 
and  only  a  slowly  varying  function  of  t: 

Z8(r)  -  Zj(t>  tt  +  Z*(t)  H  (3. 14) 

2S  will  be  discretized  over  the  coordinate  t  but  will  be 
assumed  to  be  approximately  constant  over  any  particular  segment. 
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In  applying  the  moment-method  for  the  IBC  case,  let  A  be 
expanded  similar  to  J : 


Mt(t,4)  -  ?  I  Mjk  fk(t)  e*n* 

n-1  k-1  K 

(3.15) 

M*(t,f)  -  I  I  M*k  fk(t)  e^n* 

n-1  k-1  nk  * 


Using  equations  2.9  and  3.12,  the  elements  of  [Mn]  can  be 
expressed  in  terms  of  the  [In]  as: 


<zJ)k 

(3.16) 

MJk  -  "  <2»>k  Ink 


Substitute  the  expansions  for  3  and  A,  equations  3.8  and 
3.15,  into  each  IBC  integral  equation  and  take  the  appropriate 
inner  products  with  the  testing  functions.  For  the  nttl  Fourier 
mode,  Eqs.  3.5  and  3.6  can  be  expressed  in  matrix  form 
respectively  as: 
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(3.17) 


‘V  -  ‘V  IV  +  7  I")  (V  tMnl 

1%)  -  'V  'V  +  {  ft  IV  tM„]  (3.18) 

A 

The  matrix  [N]  represents  the  n  x  vector  operation  and  can  be 
expressed  as : 


IN]  - 


I 


0 


(3.19) 


where  the  K  x  K  submatrix  I  is  the  identity  matrix.  It  is 
desired  to  convert  Eqs.  3.17  and  3.18  into  the  following  forms 
respectively: 

[Vn3  "  [ZZn]  [In]  (3.20) 

[Un]  -  [YYn]  [In]  (3.21) 

where  [ZZn]  and  [ YY^]  can  be  thought  of  as  the  generalized, 
normalized  impedance  and  admittance  matrices  for  the  IBC 
equations.  In  each  of  Eqs.  3.17  and  3.18,  there  is  a  term  of  the 
form  [N]  [Ar]  [Mn]  which  must  be  converted  into  the  form 
[Bn][In].  This  is  accomplished  by  the  following  steps. 

A 

Applying  the  n  x  operation: 
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where  the  elements  of  [Bn]  are  given  by:  (i-row,  j=column) 


<BnC>lj  '  <hi  <^>J 

<B”>u  -  -<AJc>ij  <22>j 

(3.25) 

Oij  -  -<An*)y  <Zs>j 

«'n'h,  -  <Olj  <Z2>J 


The  matrix  elements  appropriate  to  each  IBC  equation  can  now  be 
written  down  by  inspection  (in  terms  of  the  matrices  pertaining 
to  the  perfect  conductor  equations).  For  the  nc^  Fourier  mode, 
the  matrix  representations  of  Eqs.  3.5,  3.6  and  3.7  are: 


EFIE  (3.26) 


MFIE  0,27) 


u!l+aV!; 

n  n 

u*+«vj 

n  n 


YY“+oZZtt  YY^+oZZ^ 
n  n  n  n 

YY^+aZZ**  YY**+aZZ** 
n  n  n  n 


CFIE  (3.28) 
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where  the  elements  of  [ZZn]  are  given  by:  (i*row,  j»column) 


<22“>tj  -  <2“>ij  +  <2s>j 

<22S*>ij  ■  <2S*>u  -  <zi>j 

<22f>ij  *  <2r>u  -  <2s>j 

<22">lj  =  <z**)tJ  ♦  (Zj)j 

The  matrix  elements  for  [YYn]  are  the  same  as  those  for  [ZZn] 
but  with  the  roles  of  [Zn]  and  [Yn]  on  the  right  hand  side  of 
3.19  interchanged. 

As  in  the  perfect  conductor  case,  once  each  matrix  equation 
is  solved  for  the  [In],  then  that  Fourier  mode  contribution  to 
the  scattered  field  can  be  calculated.  In  the  IBC  case,  however, 
both  3  and  contribute  to  the  scattered  field.  The  next  section 
will  consider  how  to  generalize  the  far  field  calculation 
algorithms  of  Mautz  and  Harrington  to  include  the  contribution  of 

ft. 

3. A  Generalization  of  Far  Field  Calculations 
In  general,  the  scattered  electric  field  due  to  equivalent 
electric  and  magnetic  sources  (3, ft)  on  a  surface  can  be  expressed 


un  ;ij 
<Yn*>lj 


<YnC>lj 


(3.29) 


2S-  -jwy  A-VV-Vx? 


(3.30) 
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where 


/  3(r')  *  ds '  (electric  vector  potential) 

s 

-  /  -r—  [v'*5(r')]  *  ds '  (scalar  potential) 

S  J  u> 

/  fi(r')  *  ds '  (magnetic  vector  potential) 
s 


At  a  far  field  point,  defined  by  the  spherical  coordinates 

AAA 

(r,8,+)  and  unit  vectors  (r,0,$),  Eq.  3.30  can  be  approximated 


fiS(?) 


Ee  0  +  K  ♦ 


where : 


E®  *  -jkn  /_[3(r')*0  ft(r')»$]  e"jkr*r  ds 


»nr  J  s 


(3.31) 


e;  -  -JkT1  U 


Mautz  and  Harrington  define  measurement  matrices,  [Rn],  having 
the  elements: 

RjJ5  *  k  /  dt'  p’f^t')  jl*  d^'(J.q)  ej[_k  r*r'+n(*'-*)  (3.32) 


where  p'  is  the  cylindrical  coordinate  p  on  the  surface  (a 


function  of  t')  and  where: 
p  ■  t  or  t 
q  ■  0  or  t 


surface  coordinate 
spherical  coordinate 


Using  3.32  and  the  expansions  for  3  and  A,  the  far  scattered 
field  components  can  be  written  as: 


-jn 


e-JVcr 

4irr 


+  £  <1  "ni  +  {  "nil  eJn* 


(3.33) 


.“Jkr 


J  -  *BT  l  l  'Rnl  lnl+ 


1  Rt0  mC  -  1  R*0  M*  1 
n  ni  Mni  n  Kni  MniJ  e 


(3.34) 


Using  equation  3.16  for  the  elements  of  [Mn],  one  obtains: 


,-Jkr 


E«  -  -J"  W  l  |  lRnl  4+  SKi  ^l’  * 


jn* 


(3.35) 


»-Jkr 


e;  ■  -j"  v  i  \  [sSi  ini+  s*nt  *Sii eJn* 
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where 


Since  the  moment-method  programs  of  Mautz  and  Harrington 
calculate  the  elements  of  {Rn],  the  [Sn]  elements  can  easily  be 
calculated  and  the  remainder  of  the  far  field  calculation  is 
essentially  the  same  as  described  in  Ref.  12. 
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4.0  NUMERICAL  RESULTS 


4.1  Introduction  and  Definitions 

Presented  in  this  section  are  numerical  solutions  to  the  IBC 
integral  equations  for  scattering  from  impedance  spheres.  The 
primary  objective  is  to  contrast  the  EFIE,  MFIE  and  CFIE 
formulations  in  terms  of  solution  accuracy  and  stability.  A 
second  objective  is  to  determine  the  sensitivity  of  the  CFIE  to 
parameters  such  as  integral  equation  weighting  coefficient  (a) 
and  moment-method  segmentation  (SEG)  for  different  types  of 
impedance  surfaces  (i.e.  perfectly  conducting,  lossy,  reactive, 
etc.) . 

The  sphere  surfaces  are  assumed  to  be  characterized  by  a 
scalar,  constant  value  of  normalized  surface  impedance  (Zs). 
Specific  values  of  Z8  are  chosen  to  represent  various  types  of 
impedance  surfaces  that  might  be  encountered  in  scattering 
problems.1"7*9  For  example,  a  pure  real  value  of  Z8  could 
represent  the  surface  impedance  of  a  Salisbury  screen  type  of 
configuration.17  The  special  case  of  Zs**1  .0  (ideal  Salisbury 
screen)  is  of  interest  since  the  flat  plate  reflection 
coefficient  at  normal  incidence  is  zero.  A  pure  imaginary  value 
of  Z8  could  represent  the  surface  impedance  of  a  rough 
perfectly  conducting  surface  or  a  perfect  conductor  coated  with 
lossless  dielectric  layers.  Reactive  values  of  Z8  can  allow 
surface  wave  resonances  to  exist  where  the  scattering  cross 


section  might  change  by  20  to  30  dB  for  a  very  small  change  in 
frequency.  Bodies  with  finite  conductivity  or  perfect  conductors 
coated  with  lossy  layers  can  be  represented  by  values  of  Z8 
having  both  real  and  imaginary  parts.  For  example  Z8~0.1+j0.1 
could  represent  the  surface  impedance  of  a  homogeneous,  lossy 
conductor  or  a  perfect  conductor  coated  with  a  thin  lossy 
dielectric  layer.  Zs“1 .0+jl .0  might  represent  a  conductor 
coated  with  a  thick  lossy  layer. 

The  numerical  results  are  presented  in  three  parts. 

Section  4.2  contains  examples  of  scattering  from  a  perfectly 
conducting  sphere  (Z8=0.0)  for  ka  values  in  the  vicinity  of  the 
first  theoretical  sphere  interior  resonance  (ka“2.744).  These 
results,  which  are  similar  to  those  of  previous 
researchers7 > 1 1 » 12 ,  review  present  understanding  of  the 
effects  of  spurious  resonances  upon  the  stability  and  accuracy  of 
moment-method  solutions  to  the  perfect  conductor  EFIE,  MFIE  and 
CFIE  formulations.  Also,  the  results  demonstrate  general 
properties  of  the  convergence  of  moment-method  solutions  as  a 
function  of  body  segmentation  and  review  techniques  for 
determining  the  best  choice  of  a  in  the  CFIE.  In  Section  4.3, 
examples  pertaining  to  scattering  from  impedance  spheres  are 
presented  which  demonstrate  the  influence  of  non-zero  values  of 
Z8  upon  the  convergence  of  solutions  to  the  CFIE.  Also,  the 
problem  of  determining  the  best  choice  of  a  for  non-zero  values 


28 


r 


of  Z8  Is  considered  in  this  section.  Finally,  Section  4.4 
contains  numerical  examples  which  directly  compare  the  stability 
and  accuracy  of  scattering  solutions  of  the  EFIE,  MFIE  and  the 
CF1E  for  non-zero  values  of  Z8  at  ka  values  near  theoretical 
Interior  sphere  resonances.  The  examples  include  backscatter  vs 
ka  near  the  first  sphere  resonance  for  Z8"0.1+j0.1  and 
Z8-1 .0,  backscatter  vs.  Re(Zs)  or  Im(Z8)  at  ka-2.75  and 
bistatic  scattering  from  a  sphere  having  Z8“1.0  at  ka  values 
near  each  of  the  first  four  sphere  resonances. 

For  purposes  of  definition,  the  term  "stability"  of  a 
moment-method  solution  refers  to  the  potential  for  numerical 
errors  in  the  matrix  equation  which  approximates  the  integral 
equation.  A  quantitative  measure  of  solution  stability  which  is 
used  in  all  examples  is  the  matrix  condition  number  (described  in 
the  next  few  paragraphs).  The  term  "accuracy"  refers  to  the 
final  solution  (cross-section)  accuracy.  The  moment-method 
solution  accuracy  for  these  sphere  examples  is  determined  by 
comparison  with  the  Mie  series  solution  (computed  for  each 
example  using  an  independent  computer  program  based  on  Ref.  14). 

As  stated  above,  a  matrix  condition  number  is  used  as  a 
measure  of  the  stability  of  the  moment-method  solutions. 

Consider  the  matrix  equation  Ax-y  where  A  is  a  square  matrix  and 
x  and  y  are  column  matrices  (or  vectors).  When  x  is  unknown  and 
y  is  known,  the  solution  is  obtained  by  calculating  A"^y.  An 
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error  in  the  calculation  of  y  of  6y  will  produce  an  error  in  x  of 
fix.  When  6y  is  in  the  direction  of  an  eigenvector  of  A,  it 
follows  from  Ax=y  that  fix  will  be  in  the  direction  of  the  same 
eigenvector.  From  the  definition  of  matrix  and  vector  norms,  the 
relative  uncertainty  in  x  can  be  related  to  that  of  y  by 

1  IQ 

(assuming  no  uncertainty  in  the  calculation  of  A  or  A”1): 


LI 

II 


<  cond  (A) 


(4.1) 


where  cond(A),  or  condition  number  of  A,  is  defined  as: 


cond  (A)  =  A  •  A“1 


(4.2) 


and  where  | |  | |  denotes  matrix  norm  when  applied  to  A  and  vector 

norm  when  applied  to  x  or  y.  When  A  is  real  and  symmetric  and 
when  the  geometric  vector  norm  is  used  (length  of  the  vector), 
cond (A)  is  the  ratio  of  the  maximum  to  the  minimum  eigenvalues  of 
A.18  If  A  is  near  singular  (i.e.  minimum  eigenvalue  near  zero) 
cond (A)  becomes  very  large.  In  such  a  case,  an  error  fix  in  the 
direction  of  the  eigenvector  associated  with  the  minimum 
eigenvalue  may  be  large  and  perhaps  may  dominate  the  numerical 
solution  for  x.  Thus,  cond (A)  can  be  interpreted  as  an 
amplification  factor  for  numerical  errors  in  solving  a  matrix 
equation.  For  example,  a  condition  number  of  1000  means  that  the 
expected  relative  uncertainty  in  the  solution  x  is  1000  times 
that  of  y. 
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To  connect  the  above  discussion  with  the  impedance  body 
scattering  formulations,  recall  that  the  calculation  of  cross 
section  described  in  Section  3  is  a  two-step  process.  First  the 
integral  equation  (EFIE,  MFIE  or  CFIE)  is  approximated  by  a 
matrix  equation  to  obtain  a  discrete  representation  of  the 
equivalent  currents  on  the  body  surface.  Then  the  scattered 
field  due  to  these  equivalent  currents  is  computed  by  a  closed 
surface  integration.  The  condition  numbers  of  the  matrix 
approximation  to  the  IBC  integral  equations  characterize  the 
stability  of  the  first  half  of  the  calculation  of  cross-section. 
That  is,  the  condition  numbers  reflect  the  probable  errors  in  the 
equivalent  currents  -  not  necessarily  the  same  as  the  relative 
errors  in  cross  section.  For  the  EFIE  and  MFIE  formulations, 
when  at  a  frequency  near  a  spurious  resonance,  the  integral 
operators  are  near  singular  and  matrix  approximations  are 
ill-conditioned.  By  analogy  with  the  example  of  the  real, 
symmetric  matrix,  the  probable  errors  in  the  equivalent  currents 
are  the  equivalent  currents  of  the  interior  resonance  problem. 

For  the  results  in  this  report,  cond(A)  is  determined  using 
the  infinity  norm  definition  -  which  is  much  easier  to  calculate 
than  the  geometric  norm  but  has  a  more  difficult  geometrical 
interpretation.18  To  accomplish  this  calculation,  the  matrix 
decomposition  and  solution  routines  of  Ref.  13  are  replaced  by 
those  of  Ref.  19.  Also,  for  all  examples  axial  incidence  to  the 
spheres  are  assumed  (without  loss  of  generality).  In  this  case 
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only  the  M-1  Fourier  mode  of  the  moment -method  solution  is 
non-zero  and  solution  stability  is  characterized  by  a  single 
condition  number. 

For  purposes  of  quantifying  solution  accuracy  in  some 
examples,  a  mean  percentage  error  (ME)  in  the  bistatic  cross 
section  is  defined  as  follows: 


ME 


1  N 
n=1 


’cfie 


<#n>  - 


3mie 


<#n> 


0mie^6n^ 


x  100% 


In  the  above  equation,  8n  is  a  bistatic  angle  and  N  is  the 

total  number  of  bistatic  angles  used.  For  all  examples,  N=362 

(181  E-plane  and  181  H-plane  angles).  acfie  represents  the 

bistatic  cross  section  of  the  CFIE  equation  and  omie  is  the  Mie 

series  solution.  Since  the  above  definition  blows  up  when  omie 

2  -5 

becomes  very  small,  for  any  value  of  o/ira  <.10  ,  a  value  of 

10~5  is  used.  This  modification  is  only  required  for  the 
examples  of  Zs“1 .0  at  small  bistatic  angles. 

As  an  aid  to  interpreting  the  ME,  Fig.  4.1  shows  four 
examples  E-plane  bistatic  scattering  patterns  for  an  impedance 
sphere  with  Zs*1.0+j1.0  and  ka=2.75.  Different  values  of  a 
result  in  different  degrees  of  solution  accuracy.  These  bistatic 
scattering  patterns  show  that  for  a  given  ME,  larger  percentage 
errors  occur  near  backscattering  while  smaller  percentage  errors 
occur  near  forward  scattering. 
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Fig«  4.1.  Example  E-plane  bistatlc  scattering  patterns 
Illustrating  different  degrees  of  solution  accuracy  and  the 
corresponding  values  of  mean  error  (ME) . 


4.2  Perfectly  Conducting  Sphere  Examples 

This  section  contains  numerical  examples  applicable  to 
scattering  from  perfectly  conducting  spheres  (Z8«0.0).  The 
purposes  of  first  presenting  these  examples  are  to:  (1)  provide 
a  reference  for  comparing  the  scattering  examples  in  which  Zs 
is  non-zero,  (2)  illustrate  general  criteria  concerning 
convergence  of  a  moment-method  solution,  (3)  review  current 
understanding  of  the  effects  of  spurious  resonant  solutions  on 
the  EFIE  and  MFIE,  and  (4)  review  methods  for  determining  a 
suitable  choice  of  integral  equation  weighting  coefficient  (a)  in 
the  CFIE.  Some  of  these  examples  have  been  considered  in  detail 
by  Oshiro,  et.  al.7  and  Mautz,  et.  al. 11-13 

Figure  4.2  illustrates  the  effect  of  increasing 
moment-method  segmentation  (SEG)  upon  the  accuracy  and  stability 
of  solutions  to  the  EFIE,  MFIE  and  CFIE  (a*. 25)  at  ka**2.50. 
Plotted  are  the  bistatic  cross  section  mean  error  (ME)  and 
condition  number  (cond)  vs.  SEG.  An  alternate  abscissa  is 
provided  below  the  figure  to  indicate  matrix  sizes.  SEG  is 
defined  as  the  number  of  straight-line  body  profile  segments  per 
wavelength.  The  number  of  expansion  functions  per  wavelength  is 
approximately  SEG/2.  In  the  figure,  the  ordinates  that  increase 
with  increasing  SEG  are  condition  numbers  and  the  ordinates  that 
decrease  are  mean  errors,  a-. 25  is  used  for  the  CFIE  solutions 
following  guidelines  established  in  the  literature7*11-13 


Z S-0  0  KA-2 . 50 


SEGMENTS  per  wavelength 
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Fig.  4.2.  Mean  error  in  the  perfectly  conducting  sphere  bistatic 
cross  section  (decreasing  ordinates)  and  condition  number 
(increasing  ordinates)  vs.  moment-method  segmentation  or  matrix 
size  (ka  -  2.50). 


(an  example  showing  how  a  suitable  choice  of  o  can  be  determined 
is  presented  later  in  this  section). 

At  this  value  of  ka  the  EFIE,  MFIE  and  CFIE  behave  similarly 
as  SEG  is  increased  -  that  is  the  accuracy  of  iach  equation 
improves  and  the  condition  numbers  of  each  equation  increase.  The 
improvement  in  accuracy  is  due  to  the  finer  moment-method 
segmentation  allowing  a  better  approximation  to  the  equivalent 
currents  while  the  increase  in  condition  numbers  is  due  to  the 
larger  matrix  sizes  having  smaller  minimum  eigenvalues.  If  the 
segmentation  is  increased  further  and  further,  eventually  a 
crossover  point  in  improvement  in  accuracy  would  be  reached 
beyond  which  the  accuracy  would  worsen  with  increasing  SEG.20 

The  accuracy  of  the  three  equations  as  described  by  the  ME 
is  comparable  at  any  given  SEG.  However,  the  relative 
stabilities  of  the  three  equations  at  a  given  SEG  are  different. 
For  matrix  sizes  larger  than  about  10  x  10,  the  EFIE  is 
relatively  ill-  conditioned  as  compared  to  the  MFIE  and  the  CFIE 
(a=.25).  This  difference  in  conditionality  is  a  result  of  the 
EFIE  being  a  Fredholm  integral  equation  of  the  1st  kind  while  the 
MFIE  and  CFIE  (which  is  heavily  weighted  with  the  MFIE  in  this 
example)  are  second  kind  equations.21  Briefly,  as  the  matrix 
size  increases  due  to  finer  moment-method  segmentation,  matrix 
elements  obtained  from  integrals  will  decrease  in  magnitude  since 
the  surface  area  of  integration  decreases.  In  a  1st  kind 
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equation  (EFIE)  all  elements  are  calculated  from  terms  involving 
surface  integrals,  and  hence  their  magnitudes  (and  the  magnitude 
of  the  minimum  eigenvalue)  decrease  as  the  matrix  size 
increases.  As  the  matrix  size  tends  toward  infinity,  the  matrix 
elements  of  a  1st  kind  equation  tend  toward  zero.  In  a  2nd  kind 
equation  (MFIE)  a  term  exists  which  does  not  involve 
integration.  Thus,  the  MFIE  tends  to  become  more  diagonally 
dominant  as  the  matrix  size  increases  -  with  the  diagonal 
elements  in  general  tending  toward  non-zero  values  as  the  matrix 
size  increases  toward  infinity.21 

The  relative  values  of  condition  number  imply  that  the 
errors  in  the  computed  equivalent  currents  are  probably  greatest 
in  the  EFIE  and  least  in  the  MFIE.  In  fact  for  some  matrix 
sizes,  the  equivalent  current  errors  in  the  EFIE  solutions  are 
probably  an  order  of  magnitude  larger  than  those  of  the  MFIE  or 
the  CFIE  (a=.25).  Yet  the  solution  accuracy  of  the  EFIE  as 
described  by  the  ME  is  generally  better  than  the  MFIE  at  a  given 
SEG  for  ka*2.50.  This  apparent  contradiction  has  been  explained 
by  examining  the  fields  generated  by  spurious  equivalent 
currents.11*12  Resonant  cavity  equivalent  currents  of  the  form 
(3,0)  satisfy  the  boundary  condition  nxS_“0  on  S  for  the  EFIE. 
Applying  Eq.  2.1  (with  ?i=0) ,  the  resonant  cavity  equivalent 
currents  satisfy  nx£+=*0  which  implies  nxft+=0  when  radiation 
conditions  are  met.  Thus  the  spurious  equivalent  currents  which 
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enter  scattering  solutions  to  the  perfect  conductor  EFIE  have 
zero  exterior  electromagnetic  fields.  If  calculated  exactly, 
these  spurious  currents  do  not  contribute  to  the  computation  of 
scattered  field.  In  contrast,  resonant  cavity  equivalent 
currents  of  the  MFIE  of  the  form  (3,0)  satisfy  nxft_*0  -  which 
implies  nxli+«3  and  nx£+*0.  These  resonant  currents  are 
non-physical  but  never-the-less  exist  when  enforcing  nxft_«0. 
Thus,  spurious  equivalent  currents  of  the  MFIE  have  non-zero 
exterior  fields  and  hence  contribute  directly  to  the  computed 
scattered  field.  The  net  result  is  that  the  final  solution 
accuracies  of  the  EFIE  and  the  MFIE  are  comparable  even  though 
the  computed  equivalent  currents  of  the  EFIE  contain  much  larger 
errors  than  those  of  the  MFIE  (verified  by  Mautz  through  actual 
computation  of  currents  in  the  sphere  examples). 

Figures  A. 3  a  and  b  illustrate  the  relative  accuracies  and 
stabilities  respectively  of  solutions  to  the  EFIE,  MFIE  and  CFIE 
(a“.25)  for  a  range  of  ka  values  near  the  first  theoretical 
sphere  resonance  (ka**2.744).  Plotted  in  Fig.  4.3a  is  normalized 
backscatter  (o/ira^)  from  a  perfectly  conducting  sphere  vs.  ka 
as  predicted  by  each  of  the  integral  equations  and  the  Mie 
series.  The  matrix  size  is  26  x  26  for  all  points  which 
corresponds  to  a  SEG  of  approximately  20/x.  Figure  4.3a 
demonstrates  the  regions  of  inaccuracy  in  the  EFIE  and  MFIE  and 
the  potential  errors  in  magnitude  in  computing  backscatter  at 
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Fig.  4.3a.  Normalized  backscatter  (o/wa^)  from  a  perfectly 
conducting  sphere  vs.  normalized  wavenumber  (ka)  for  ka  values 
near  the  first  theoretical  sphere  interior  resonance  (ka  ■ 
2.744).  Segmentation  (SEG)  is  defined  as  the  number  of 
straight-line  body  profile  segments  per  wavelength. 


SEG-20A.  Plotted  in  Fig.  4.3b  are  the  condition  numbers  which 
correspond  to  the  data  of  Fig.  4.3a.  These  values  of  cond 
indicate  the  regions  of  instability  and  the  relative  errors  in 
the  equivalent  currents  of  solutions  to  the  three  integral 
equations. 

Figures  4.3a  and  b  both  indicate  that  the  EFIE  and  MFIE  are 
affected  by  (or  predict  the  location  of)  the  spurious  resonance 
at  different  ka  values.  By  virtue  of  its  inherent  better 
conditioning,  the  MFIE  predicts  the  resonance  closer  to  the 
theoretical  value.  Figures  4.3a  and  b  also  indicate  that  the  1st 
sphere  interior  resonance  affects  the  MFIE  over  a  broader  range 
of  ka  values  than  the  EFIE  (i.e.,  the  MFIE  has  a  wider  bandwidth 
for  errors).  However,  the  largest  potential  magnitude  of  errors 
in  backscatter  from  the  sphere  appears  to  occur  for  the  EFIE  In  a 
very  narrow  range  of  ka  values.  Throughout  the  region  of  the 
interior  sphere  resonance,  the  CFIE  (a*. 25)  gives  accurate 
scattering  predictions  and  demonstrates  no  ill-conditioning  from 
the  resonance. 

Figure  4.4  shows  the  effect  of  changing  matrix  size  on  the 
scattering  predictions  of  the  EFIE  near  the  first  sphere 
resonance.  In  this  figure,  normalized  backscatter  from  a 
perfectly  conducting  sphere  vs.  ka  is  plotted  using  three 
different  matrix  sizes  in  the  EFIE  solution  -6x6  (SEG-7/A), 
12x12  (SEG-10/A)  and  26  x  26  (SEG-20/A).  The  dashed  ordinate 
is  the  Mie  series  solution.  As  matrix  size  Increases,  three 
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Fig.  4.4.  Effect  of  changing  moment-method  segmentation  on  the 
normalized  backscatter  vs.  ka  predictions  of  the  EFIE  (perfectly 
conducting  sphere). 


effects  can  be  seen  upon  the  contamination  of  the  EFIE  scattering 
solutions  by  the  spurious  resonance.  First,  the  magnitude  of  the 
potential  errors  decreases  as  matrix  size  increases.  Secondly, 
the  region  of  inaccuracy  (bandwidth  of  errors)  narrows  as 
indicated  by  the  shape  of  the  curves.  Thirdly,  the  predicted 
location  of  the  resonance  gets  closer  to  the  theoretical  ka  value 
as  the  matrix  size  increases.  The  MFIE  scattering  predictions 
(not  shown  in  the  figure)  are  affected  in  a  similar  manner  -  with 
the  bandwidths  for  errors  being  wider  than  those  of  the  EFIE  at 
the  same  matrix  size  and  the  predicted  location  of  the  resonance 
being  closer  to  the  theoretical  value  at  a  given  matrix  size. 

The  problem  of  choosing  a  suitable  value  for  the  CFIE 
weighting  coefficient  (a)  has  been  considered  by  Oshiro,  et.al.7 
and  Mautz,  et.al.11'13  for  scattering  from  perfectly  conducting 
spheres.  Using  a  bistatic  cross  section  mean  error  (very  similar 
to  the  one  used  in  this  report),  Oshiro  determined  the  best 
choice  of  a  to  be  in  the  range  of  .2  -  .3  for  a  variety  of  ka 
values.  Using  the  rms  error  in  the  equivalent  current  as  the 
criterion  (analogous  to  the  condition  numbers  of  this  report), 
Mautz  and  Harrington  determined  the  best  choice  of  a  to  be 
approximately  .2  for  ka  near  the  1st  sphere  resonance. 

Figure  4.5  illustrates  the  above  two  described  methods  by 
displaying  ME  and  cond  vs.  a  for  5  different  ka  values  near  the 
1 8 1  sphere  resonance.  Matrix  sizes  are  26  x  26  which  corresponds 
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to  SEG-20/X.  When  a-. 01,  the  CFIE  solution  can  be  considered 
near  to  the  MFIE  solution.  When  a-100.0,  the  CFIE  solution  is 
near  an  EFIE  solution.  The  values  of  ME  are  all  <  10%  - 
indicating  that  for  SEG=20/X  good  accuracy  in  the  cross  section 
can  be  obtained  using  almost  any  a  (for  perfectly  conducting 
spheres).  It  appears  that  the  range  of  .1  <  a  <  1.0  may  be  the 
best  range  to  obtain  good  accuracy  for  any  ka.  The  values  of 
cond  in  Fig.  4.5  indicate  that  there  exists  a  more  narrow  range 
of  a  values  for  obtaining  the  most  well-conditioned  matrices  and 
smallest  errors  in  the  currents  for  all  ka  values  of  Fig.  4.5. 
Clearly,  .2  <  o  <  .3  is  a  range  which  gives  well-conditioned 
matrices  for  all  ka  values  near  the  1st  sphere  resonance. 

4.3  Influence  Of  Non-Zero  Values  Of  Surface  Impedance  On  Choice 

Of  Integral  Equation  Weighting  Coefficient  And  Moment-Method 

Segmentation  For  The  Combined  Field  Integral  Equation 

The  problems  of  choosing  the  CFIE  weighting  coefficient  (a) 
and  adequate  moment-method  segmentation  (SEG)  when  calculating 
scattering  from  impedance  spheres  are  considered  in  this  section. 
The  values  of  normalized  surface  impedance  used  for  these 
numerical  examples  are  Zs”0.1+j0.1  (representing  a  thin  lossy 
dielectric  layer),  Zs=1.0+j1.0  (representing  a  thick  lossy 
layer),  Z8->±j1.0  (representing  lossless  inductive  and 
capacitive  surfaces  and  Z8*1.0  (representing  an  ideal  Salisbury 


screen).  Choosing  the  best  a  is  considered  first  for  2  ka  values 
near  the  1st  sphere  resonance  assuming  that  SEG-20/X  is  adequate 
segmentation.  Using  suitable  choices  of  a  for  the  different 
impedance  surfaces,  the  convergence  of  the  CF1E  as  a  function  of 
SEG  is  then  examined  at  ka=2.75. 

Figures  4.6a  and  b  show  respectively  the  accuracy  and 
stability  of  the  CFIE  scattering  solutions  as  a  function  of  a  at 
ka  2.50.  ME  vs.  a  is  displayed  in  Fig.  4.6a  and  cond  vs.  i  is 
displayed  in  4.6b.  Matrix  sizes  of  the  data  in  both  figures  are 
26  x  26  (SEG-22/X).  For  the  values  of  Z8  considered,  the 
solution  accuracy  as  described  by  the  ME  is  highly  variable  for  a 
given  matrix  size.  The  best  ME  that  could  be  achieved  for 
Zg^tjl.O  is  about  7.2%  -  significantly  worse  than  the  perfectly 
conducting  case  of  .062%  (Fig.  4.2)  at  the  same  SEG.  The  values 
of  ME  in  Fig.  4.6a  also  demonstrate  that  the  functional 
dependence  of  solution  accuracy  on  a  is  also  highly  variable 
depending  on  Zs.  For  example,  values  of  ME  for  Z8~±j1.0  are 
relatively  insensitive  to  a  while  values  of  ME  when  Z8~1.0  are 
very  sensitive  to  a.  For  this  latter  case,  ME  less  than  10%  is 
obtained  only  in  a  very  narrow  range  of  o  around  1 .0. 

The  condition  numbers  of  Fig.  4.6b  show  that  for  a  given 
matrix  size  and  ka,  the  dependence  of  cond  on  a  and  the  values  of 
cond  at  a  given  o  both  appear  to  be  dependent  on  Z8.  Only  the 
case  of  Z8**0  • 1^j0. 1  was  similar  to  the  perfectly  conducting 
case  -  where  small  a  implied  well-conditioned  matrices  (cond«<20) 
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and  large  a  implied  ill-conditioned  matrices  (cond  300  -  400). 

The  influence  of  Zs  on  the  dependence  of  cond  with  a  can  be 
partly  explained  by  observing  that  Zs  plays  a  similar  role  to 
that  of  when  Z8-0  (i.e.,  fi-0) .  That  is,  the  ill-conditioned 
Lg  operator  (perfect  conductor  EFIE)  is  combined  with  the 
well-conditioned  Lh  operator  (perfect  conductor  MFIE)  in  both 
the  the  EFIE  and  MFIE  (see  Eqs.  3.5  and  3.6).  These  two 
equations  are  then  combined  in  the  CFIE.  However,  when  Z8  has 
an  absolute  value  on  the  order  of  1.0,  the  mix  of  the  two  perfect 
conductor  operators  is  about  the  same  in  the  CFIE  regardless  of 
a.  Thus,  the  value  of  Zs  (or  the  type  of  impedance  surface) 
also  plays  a  role  in  the  variability  of  the  condition  number  at  a 
given  a. 

figures  4.7a  and  b  display  ME  vs.  a  and  cond  vs.  a 
respectively  for  ka=2.75.  Matrix  sizes  are  26  x  26  - 
corresponding  to  SEG=20/X.  At  this  ka,  very  close  to  the  1st 
sphere  resonance,  the  sensitivity  to  a  of  the  CFIE  solution 
accuracy  and  stability  is  noticeably  increased  over  that  at 
ka»2.50  (Figs.  4.6  a  and  b).  Also,  the  best  solution  accuracy 
that  is  achieved  for  the  pure  reactive  surface  impedance  cases 
(Zs»±j1.0)  is  noticeably  worse  at  ka-2.75  over  that  at 
ka-2.50.  For  these  pure  reactive , cases ,  the  best  ME  at  SEG~20/X 
is  only  about  10%.  To  obtain  ME  <  10%,  choice  of  a  is  critical 
for  Z8«1.0  and  Z8-±j1.0.  For  Z8-.1+j.1  and  1.0+jl.O, 
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CflE  K  A • a  7  ?  MATRIX  SIZES*  ?  6  X  2  6 


-  LEGEND  - 
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Fig.  4.7b.  Matrix  condition  numbers  (m  -  1  Fourier  mode)  for  the 
impedance  sphere  solutions  of  Fig.  4.7a  vs.  CFIE  weighting 
coefficient  (a). 
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the  choice  of  a  is  less  critical  to  obtain  ME  <  10%.  However, 
even  in  these  latter  two  cases,  the  choice  of  a  is  more  important 
than  in  the  perfectly  conducting  case  (Fig.  4.5). 

The  final  numerical  example  in  this  section  considers  the 
influence  of  non-zero  surface  impedance  on  the  convergence  of 
impedance  sphere  scattering  solutions  of  the  CFIE.  Plotted  in 
Figs.  4.8a  and  b  respectively  are  ME  vs.  SEG  and  cond  vs.  SEG  at 
ka**2.75  for  the  different  values  of  Z8.  a  is  chosen  to  obtain 
close  to  the  best  solution  accuracy  at  a  given  Z8  using 
Fig.  4.7a.  The  actual  values  of  o  used  in  Figs.  4.8a  and  b  are: 


Za _  alpha 


0.1+j0.1 

0.6 

1  .0 

1.0 

jl.o 

7.0 

-jl  .0 

0.14 

1 . 0+ j 1  . 0 

1.25 

The  convergence  of  the  CFIE  as  a  function  of  moment-method 
segmentation  is  similar  to  the  perfect  conductor  case  (Fig.  4.2) 
in  that  increasing  SEG  generally  improves  solution  accuracy.  As 
SEG  increases  from  5/X  to  1 0/ A  there  is  a  rapid  improvement  in 
solution  accuracy.  For  SEG  >  20/X  improvement  in  accuracy  is 
3 low.  Convergence  of  the  CFIE  for  non-zero  Z8  is  different 
from  the  perfectly  conducting  case  in  that  a  given  value  of  SEG 
does  not  guarantee  that  solution  accuracy  will  be  comparable  to 
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Fig.  4.8a.  Mean  error  in  Che  impedance  sphere  bistatic  cross 
section  vs.  moment-method  segmentation  at  ka  -  2.75.  Note: 

Z8  ■  ±j1.0  ordinates  nearly  overlay  each  other,  apparently  a 
consequence  of  <*(Zg  -  jl.O)  -  1/a(Za  -  —  j  1.0). 


ZS-  0  . 1 ♦ JO . 1 
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Fig.  4.8b.  Matrix  condition  numbers  (m  -  1  Fourier  mode)  of  the 
impedance  sphere  solutions  of  Fig.  4.8a  vs.  matrix  size. 


that  of  the  perfectly  case.  This  fact  could  also  have  been 
inferred  from  Figs.  4.6a  and  4.7a  for  which  all  data  points  are 
at  SEG- 20/ A. 

The  convergence  of  the  CFIE  when  Za~±j1.0  is  interesting 
in  two  respects.  First,  the  solutions  are  relatively  inaccurate 
and  the  matrices  are  relatively  ill-conditioned  when  Z8*±j1.0 
as  compared  to  the  other  cases  of  Zs.  Secondly,  for  SEG  >  20/ A 
there  is  actually  an  increase  in  mean  errors  with  increasing  SEG 
while  for  the  other  values  of  Zs  there  is  a  gradual  decrease  in 
ME.  The  increase  in  ME  with  increasing  SEG  for  the  pure  reactive 
cases  suggests  that  the  best  choice  of  a  changes  when  SEG  is 
changed.  However,  this  possibility  was  determined  not  to  be  the 
reason  for  the  convergence  behavior  for  Z8-±j1.0  exhibited  in 
Fig.  4.8a.  The  best  choice  of  a  was  determined  for  Z8»j1.0  at 
SEG=32/a  (matrix  sizes=*45  x  45)  by  generating  a  curve  similar  to 
the  one  for  SEG=20/X  in  Fig.  4.7a.  The  results  were  that  the 
best  choice  of  a  remained  at  approximately  7.0  with  the  lowest  ME 
approximately  15%.  Thus,  it  appears  that  the  cross-over  point  at 
which  there  is  no  longer  improvement  in  accuracy  with  increasing 
matrix  size  occurs  very  early  for  Z8*±j1.0  at  ka-2.75. 
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4.4  Comparison  Of  the  Electric,  Magnetic  and  Combined  Field 

Integral  Equations  For  Non-Zero  Values  Of  Surface  Impedance 

Examples  of  scattering  from  impedance  spheres  are  presented 
in  this  section  which  directly  compare  solutions  to  the  EFIE  and 
MFIE  with  those  of  the  CFIE  for  ka  values  near  spurious 
resonances.  The  primary  objective  of  this  section  is  to 
demonstrate  that  spurious  resonant  solutions  can  severely 
contaminate  the  EFIE  and  MFIE  predictions  of  exterior  scattering 
from  impedance  bodies  while  the  CFIE  is  a  formulation  which 
rejects  the  spurious  solutions. 

Backscattering  from  a  sphere  with  with  Z8»0.1+j0.1  at  ka 
values  near  the  first  sphere  resonance  is  considered  in 
Figs.  4.9a  and  b.  Normalized  backscatter  vs.  ka  as  predicted  by 
the  EFIE,  MFIE,  CFIE  (a=.6)  and  Mie  series  is  plotted  in 
Fig.  4.9a  and  condition  numbers  for  the  integral  equations  vs.  ka 
in  Fig.  4.9b.  Matrix  sizes  of  the  data  are  26  x  26  which 
corresponds  to  SEG-20/A.  The  scale  sizes  of  Figs.  4.9a  and  b  are 
the  same  as  those  of  Figs.  4.3a  and  b  to  allow  direct  comparison 
with  the  perfectly  conducting  case. 

The  presence  of  the  non-zero  surface  impedance  causes  some 
subtle  changes  in  the  accuracy  and  stability  of  the  integral 
equation  solutions  relative  to  the  perfectly  conducting  case. 
First,  the  bandwidth  for  errors  and  the  magnitude  of  errors  in 
the  EFIE  and  MFIE  cross  sections  are  larger  than  in  the  perfectly 
conducting  case.  This  is  just  the  opposite  that  one  would  expect 
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Fig.  4.9a.  Normalized  backscatter  (o/wa^)  from  a  lossy 
impedance  sphere  having  a  normalized  surface  impedance  of  0.1  + 
jO.1  vs.  normalized  wavenumber  (ka). 
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if  one  assumed  that  the  lossy  surface  impedance  would  help 
regularize  the  ill-posed  equations.  Values  of  condition  number, 
however,  are  not  greatly  different  from  the  perfectly  conducting 
case  except  that  the  EFIE  condition  numbers  are  slightly  lower.  A 
possible  explanation  for  these  changes  in  accuracy  and  stability 
is  offered  as  follows.  As  mentioned  in  the  previous  section,  the 
non-zero  value  of  Zs  adds  a  term  (1 .0/n)nxLn(ft)  to  the 
perfect  conductor  EFIE  -  changing  it  from  a  1st  kind  equation  to 
a  second  kind  equation.  Hence,  one  would  expect  that  the 
conditioning  of  the  EFIE  for  this  case  should  be  better  than  the 
perfectly  conducting  case  because  of  the  addition  of  some  well- 
conditioned  Lh  integral  operator  to  the  equation.  However,  the 
non-zero  Zs  also  changes  the  spurious  equivalent  currents  from 
the  form  (3,0)  for  the  perfectly  conducting  case  to  the  combined 
source  form  (3 ,fl=-nxZg • n3) . 1 0  These  latter  currents,  which  are 
still  equivalent  currents  of  the  perfect  conductor  cavity 
problem,  always  have  non-zero  exterior  fields.  Thus,  the 
spurious  currents  of  the  EFIE  radiate  an  exterior  field  when  Zs 
is  non-zero  and  the  accuracy  in  the  cross  section  would  be 
expected  to  be  worse  than  in  the  perfect  conductor  EFIE  case.  The 
increased  errors  in  the  MFIE  cross  sections  may  be  due  to  the 
combined  source  spurious  equivalent  currents  producing  larger 
exterior  fields  than  in  the  perfectly  conducting  MFIE  case. 
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Effects  of  a  pure  real  surface  impedance  on  the  accuracy  and 
stability  of  the  three  integral  equations  is  shown  in  Fig.  4.10. 
Normalized  backscatter  (o/ira^)  and  condition  number  (cond)  are 
plotted  as  functions  of  Re(Zs)  at  ka-2.75.  Im(Zs)«=0.0  and 

matrix  sizes  of  26  x  26  are  used  (SEG=20/X).  Since  o=1.0  is  the 
best  choice  for  Zs“1.0  (as  determined  in  the  previous  section), 
this  value  is  used  for  the  CFIE  solution  at  all  values  of  Z8  in 
Fig.  4.10.  At  this  value  of  ka  (very  close  to  the  first 
theoretical  sphere  resonance),  the  EFIE  and  MFIE  are 
characterized  by  ill-conditioned  matrices  and  poor  accuracy  in 
predicting  cross  section.  In  contrast,  the  CFIE  (o=1 .0)  is 
relatively  well-conditioned  and  predicts  the  cross  section  very 
close  to  the  Mie  series  at  all  values  of  Re(Zs).  As  in  the 
previous  case,  the  lossy  surface  impedance  did  not  help  to 
regularize  the  ill-posed  integral  equation  formulations. 

The  case  of  Zs~1.0  appears  to  be  a  stressing  case  in  the 
calculation  of  scattering  from  impedance  BORs  since  the 
theoretical  cross  section  for  axial  backscattering  is  zero. 
Figures  4.11a  and  b  consider  this  case  for  a  sphere  by  displaying 
normalized  backscatter  (o/ua^)  vs.  ka  and  cond  vs.  ka 
respectively.  In  Fig.  4.11a,  there  is  a  scale  change  on  the 
vertical  axis  -  the  upper  scale  for  the  EFIE-MFIE  ordinates  and 
the  lower  scale  for  the  CFIE  (a-1.0)  ordinates.  Since  the 
backscattering  predictions  of  the  EFIE  and  the  MFIE  are 
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NORMALIZED  backscattcr  condition  NUMBER 


Fig.  4.10.  Normalized  backscatter  (lower  set  of  ordinates)  from 
a  lossy  impedance  sphere  and  associated  matrix  condition  numbers 
(upper  set  of  ordinates)  vs.  real  part  of  normalized  surface 
impedance  (imaginary  part  of  Zs  -  0.0,  ka  -  2.75).  Note:  the 
CFIE  and  Mie  series  ordinates  nearly  overlay  each  other. 


normalized  BACKSCATTER 
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Fig.  4.11a.  Normalized  backscatter  vs.  ka  for  a  lossy  impedance 
sphere  having  normalized  surface  impedance  of  1.0.  Note  the 
scale  change  on  the  vertical  axis.  Also,  the  increased 
jaggedness  of  the  CFIE  ordinates  near  the  center  of  the  figure  i 
due  to  increased  samples  in  that  region. 


63 


degenerate  for  Zs“1.0,  there  is  a  common  ordinate  for  both 
equations.  In  the  figures,  two  different  matrix  sizes  are  used 
which  illustrate  the  effect  of  a  change  in  moment-method 
segmentation  upon  the  influence  of  the  spurious  resonance  in  the 
impedance  sphere  scattering  solutions  with  Zs=1.0.  Case  1 
corresponds  to  matrix  sizes  of  26  x  26  (SEG-20/X)  and  Case  2 
corresponds  to  matrix  sizes  of  54  x  54  (SEG*>40/X). 

Considering  that  the  exact  scattering  solution  for  all  ka  is 
zero,  the  EFIE-MFIE  predictions  are  very  poor  for  all  ka  values 
considered  in  Fig.  4.11a.  The  cross  section  values  indicated  by 
the  EFIE-MFIE  ordinates  show  only  the  contribution  of  the 
spurious  equivalent  currents  to  the  computed  scattered  field.  The 
CFIE  (a=1.0)  predictions,  which  appear  to  be  computational  noise, 
are  a  much  better  approximation  to  the  theoretical  cross 
section.  The  increase  in  matrix  size  upon  the  EFIE-MFIE 
solutions  has  similar  effects  as  those  for  the  perfectly 
conducting  case  (Fig.  4.4).  Increasing  the  matrix  size  causes 
the  bandwidth  for  errors  in  cross  section  to  narrow  slightly,  the 
scattering  predictions  to  be  generally  better  at  a  given  ka  and 
the  predicted  location  of  the  resonance  to  be  shifted  closer  to 
the  theoretical  location  (ka=2.744).  However,  increasing  the 
matrix  size  has  little  effect  on  the  CFIE  (a-1.0)  other  than  an 
increase  in  the  condition  numbers  as  shown  in  Fig.  4.11b.  A 
further  improvement  in  the  CFIE  (a-1.0)  cross  section  predictions 
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would  probably  require  Increasing  computational  accuracy  (which 
is  32  bit  in  all  numerical  examples  of  this  report). 

Figures  4.12a,  b  and  c  display  E-plane  bistatic  scattering 
from  an  impedance  sphere  having  Z8»1.0  at  ka  values  near  the 
first,  second  and  third  perfect  conductor  sphere  resonances 
respectively.  In  each  case,  matrix  size  is  adjusted  to  obtain 
SEG-20/X.  Bistatic  scattering  angles  of  0  and  180  degrees 
correspond  to  backscattering  and  forward  scattering 
respectively.  The  effect  of  the  spurious  solutions  upon  the  EFIE 
and  MFIE  appears  to  be  similar  for  each  of  these  resonances  - 
with  the  predicted  location  of  the  resonance  being  slightly 
higher  than  the  theoretical  location  in  each  case  (the  ka  values 
of  each  are  very  close  to  the  peak  of  the  resonance).  In  all 
three  cases,  the  EFIE  and  MFIE  scattering  predictions  are  very 
poor  for  all  bistatic  angles  except  those  near  forward 
scattering.  The  CFIE  (o=1.0)  predictions  on  the  other  hand,  are 
excellent  for  all  bistatic  aspect  angles  -  the  curves  nearly 
overlaying  the  Mie  series  solutions  in  all  three  cases.  Although 
o=1.0  was  determined  in  the  previous  section  to  be  the  best  value 
only  at  ka=2.75,  the  excellent  agreement  of  the  CFIE  (o=1.0) 
predictions  with  the  Mie  series  at  the  other  ka  values  indicates 
that  a=1.0  may  be  the  best  value  at  all  frequencies  when  the 
normalized  surface  impedance  is  1.0. 
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Fij».  4.12a.  Normalized  E-plane  bistatic  scattering  from  a  lossy 
impedance  sphere  (Zg  -  1.0,  ka  -  2,758).  For  the  H-p lane, 
interchange  the  EFIE  and  MFIE  predictions  (the  CFIE  and  Mie 
series  predictions  for  the  H-plane  are  the  same  as  those  for  the 
E-plane  when  Z8  ■  1.0). 
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Fig.  4.12c.  Normalized  E-plane  bistatic  scattering  from  a  lossy 
impedance  sphere  (Za  «  1.0,  ka  *  4.51). 


Figures  4.13  and  4.14  consider  scattering  from  impedance 
spheres  characterized  by  a  pure  reactive  surface  impedance. 

Figure  4.13  displays  normalized  backscatter  (o/va^)  vs. 

Im(Zs)  at  ka=2.75  with  Re(Zs)**0.  Four  different  integral 
equation  solutions  are  compared  to  the  Mie  series:  EFIE,  MFIE, 
CFIE  (a*1 . 0)  and  CF1E  (variable  a).  The  last  case  is  presented 
to  demonstrate  that  an  improvement  (albeit  slight)  in  accuracy 
can  be  attained  by  varying  o  as  Z8  changes.  Figure  4.14 
displays  condition  number  vs.  Im(Z8)  for  each  of  the  integral 
equation  solutions  of  Fig.  4.13.  Matrix  sizes  are  26  x  26  - 
corresponding  to  SEG-20/X. 

The  accuracies  of  the  EFIE  and  the  MFIE  in  predicting 
backscatter  at  ka-2.75  are  very  poor  when  Z8  is  pure 
imaginary.  In  contrast,  both  CFIE  solutions  have  good  agreement 
with  the  Mie  series.  The  CFIE  (variable  a)  case  represents  a 
slight  improvement  over  the  constant  a  case  -  although  there  are 
still  large  errors  in  cross  section  in  the  regions  where  cross 
section  changes  greatly  with  a  slight  change  in  Im(Z8).  Values 
of  condition  number  displayed  in  Fig.  4.14  Indicate  that  the  CFIE 
solutions  are  generally  better  conditioned  than  the  EFIE  and  MFIE 
solutions.  The  stability  of  the  integral  equation  solutions  for 
large  inductive  Z8  appears  to  be  strikingly  different  than  that 
for  large  capacitive  values  of  Z8.  That  is,  the  impedance 
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Fig.  4.13.  Normalized  backscatter  from  a  reactive  impedance 
sphere  vs.  imaginary  part  of  normalized  surface  impedance  (real 
part  of  Zg  -  0.0,  ka  -  2.75).  For  the  CFIE  (variable  o)  case: 
-5.0  <  Im  (Zs)  <  -1.4,  o  -  1.0;  -1.4  <  1m  (Zs)  <  -0.6,  a  - 
.12;  -0.6  <  Im  (Z8)  <  0.6,  a  -  1.0;  0.6  <  Im  (Z8)  <  1.7,  a  - 
8.0;  and  1.7  <  Im  (Z8)  <  5.0,  a  -  1.0. 
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spheres  appear  to  be  characterized  by  resonances  for  large 
inductive  Zs  -  possibly  exterior  surface  wave  resonances.  This 
behavior  is  similar  to  that  noted  by  Wait  in  Mie  series  solutions 
of  scattering  from  impedance  spheres.14 
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5.0  CONCLUSIONS 

This  report  has  investigated  numerical  solutions  to  integral 
equation  formulations  for  scattering  from  bodies  described  by  a 
Leontovich  impedance  boundary  condition  (IBC).  A  method  of 
constructing  moment-method  solutions  to  three  different  IBC 
integral  equations  was  described.  Two  of  these  equations  -  the 
principal  ones  for  which  numerical  solutions  have  been  obtained  - 
are  the  electric  and  magnetic  field  integral  equations  (EFIE  and 
MFIE).  These  two  equations  are  ill-posed  because  they 
theoretically  allow  spurious  solutions  at  frequencies  at  interior 
perfect  conductor  cavity  resonances.  The  IBC  combined  field 
integral  equation  (CFIE)  -  which  has  received  little  attention  in 
terms  of  actual  numerical  solutions  -  is  well-posed  because  it 
rejects  the  spurious  resonances.  Numerical  solutions  to  the 
EFIE,  MFIE  and  CFIE  for  scattering  from  Impedance  spheres  were 
presented  in  this  report  which  considered  different  types  of 
impedance  surfaces  (i.e.,  perfectly  conducting,  lossy,  reactive, 
etc. ) 

The  numerical  results  demonstrated  that  scattering  solutions 
to  the  EFIE  and  MFIE  can  be  severely  contaminated  by  spurious 
solutions  when  at  a  frequency  near  an  interior  resonance.  This 
contamination  occurs  regardless  of  whether  the  surface  impedance 
is  lossy  or  lossless  and  no  regularization  of  the  ill-posedness 
of  these  equations  occurs  when  the  surface  Impedance  is  lossy.  In 
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contrast  to  the  MFIE  and  the  EFIE,  solutions  to  the  CFIE  are  not 
affected  by  spurious  resonant  solutions  when  a  suitable  choice  of 
integral  equation  weighting  coefficient  (a)  can  be  determined. 

The  CFIE  solutions  are  consistently  more  stable  and  accurate  than 
those  of  the  EFIE  or  the  MFIE  when  at  a  frequency  near  an 
interior  resonance.  For  some  values  of  Z8  (for  example 
Zs“1 .0)  the  solution  accuracy  of  the  CFIE  can  be  much  better 
than  that  of  the  EFIE  or  MFIE  even  when  at  a  frequency  that 
appears  to  be  far  from  an  interior  resonance.  Thus,  the 
numerical  examples  clearly  demonstrated  that  the  CFIE  is  in 
practice  (as  in  theory)  a  superior  IBC  integral  equation 
formulation  to  either  the  EFIE  or  the  MFIE. 

The  numerical  examples  also  illustrated  some  important 
differences  when  the  scatterer  surface  has  a  non-zero  surface 
impedance  than  when  the  scatterer  is  a  perfect  conductor  (this 
latter  case  being  the  one  upon  which  current  understanding  of  the 
effects  of  spurious  solutions  on  an  integral  equation  formulation 
is  based).  First,  when  Z8  is  non-zero,  relative  errors  in  the 
computed  EFIE  and  MFIE  cross  sections  are  generally  larger  and 
the  bandwidths  for  errors  are  generally  wider  than  those  for  the 
perfect  conductor  EFIE  and  MFIE.  This  appears  to  be  due  to  the 
fact  that  when  Z8  is  non-zero,  the  spurious  equivalent  currents 
take  the  combined  source  form  (3 ,ft--nx2s»n3)  rather  than  the 
form  (3,0).  The  combined  source  equivalent  currents  (which  are 
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still  equivalent  currents  of  the  perfect  conductor  interior 
resonance  problem)  radiate  an  exterior  electromagnetic  field 
which  depends  in  magnitude  upon  the  value  of  Zs.  Thus,  the  use 
of  a  well-posed  integral  equation  formulation  takes  on  greater 
importance  when  the  scatterer  surface  impedance  is  non-zero  than 
for  the  perfect  conductor  case.  Secondly,  when  Zs  is  non-zero, 
choice  of  moment-  method  segmentation  (SEG)  and  integral  equation 
weighting  coefficient  (a)  for  the  CFIE  appears  to  be  more 
complicated  to  attain  a  specified  cross  section  accuracy  than  the 
perfect  conductor  case.  Convergence  as  a  function  of  SEG  appears 
to  be  very  poor  for  certain  values  of  reactive  surface  impedance 
(low  loss)  at  certain  frequencies.  In  particular,  when  the 
frequency  is  near  an  exterior  surface  wave  resonance  (where  there 
might  be  a  very  large  cross  section  or  one  that  changes  rapidly 
with  frequency),  the. errors  may  be  large  regardless  of 
segmentation.  When  varying  a  at  a  given  SEG,  the  best  accuracy  of 
the  CFIE  occurs  at  different  values  of  a  for  different  values  of 
Zs.  Thus,  the  established  perfect  conductor  criteria  for 
choosing  a  is  not  adequate  for  all  Z8. 


APPENDIX  A 


The  Fortran  IV  computer  program  listed  in  this  appendix  is 
essentially  Mautz  and  Harrington's  main  program13  for  the  perfect 
conductor  case  extended  to  the  impedance  boundary  condition 
case.  Several  variables  (ZS,  ZST,  ZSP,  YY,  ZZ,  TM,  PM)  and  one 
subroutine  (MATB)  were  added  to  accomplish  the  required 
modifications.  Also,  many  comments  were  added  to  aid  in 
understanding  the  logic  of  the  program  and  all  the  main 
variables.  Some  minor  changes  to  Mautz  and  Harrington's  program 
were  that  the  inputs  were  changed  to  list-directed  read  and  all 
format  statements  for  outputs  were  moved  to  the  end  of  the 
program.  The  output  was  modified  to  include  printing  of  the  two 
components  of  magnetic  surface  current  density  (sampled  at  the 
same  body  locations  as  the  electric  surface  current  density). 

The  listing  does  not  include  any  routines  from  the  Mautz  and 
Harrington  program  that  are  unchanged. 

Solution  is  obtained  for  the  following  conditions:  (1) 
single  excitation  frequency,  (2)  one  incident  angle  and  one 
scattering  angle  (not  necessarily  the  same),  and  (3)  one  value  of 
normalized  surface  impedance  on  the  BOR.  The  actual  program  used 
for  generating  the  results  of  Chapter  4  was  a  modification  of  the 
one  listed  in  this  appendix  which  could  calculate  scattering  at 
many  aspect  angles,  frequencies,  etc.  (by  appropriate  insertion 
of  loops ) . 
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PROGRAM  IBC  -  THIS  PROGRAM  SOLVES  FOR  THE  ELECTROMAGNETIC 
SCATTERING  FROM  A  BODY  OF  REVOLUTION  HAVING  A  SURFACE 
CHARACTERIZED  BY  A  LEONTOVICH  IMPEDANCE  BOUNDARY  CONDITION 
(IBC) »  A  MOMENT- METHOD  SOLUTION  IS  OBTAINED  FROM  THE 
ELECTRIC  FIELD  INTEGRAL  EQUATION  (EFIB) ,  THE  MAGNETIC 
FIELD  INTEGRAL  EQUATION  (MFIZ)  AND  THE  COMBINED  FIBLD 
INTEGRAL  EQUATION  (CFIE)  .  THE  INTERNAL  CODING  OF  THE 
PROGRAM  ASSUMES  A  NORMALIZED,  DIAGONAL,  DYADIC  SURFACE 
IMPEDANCE  HHICH  CAN  HAVE  ARBITRARY  VARIATIONS  ALONG  THE 
LONGITUDINAL  SURFACE  COORDINATE  OF  THE  BOR.  HOWEVER,  THIS 
VERSION  READS  IN  A  SCALAR,  CONSTANT  VALUE  OF  NORMALIZED 
SURFACE  IMPEDANCE* 

THIS  PROGRAM  IS  A  GENERALIZATION  MAUTZ  AND  HARRINGTON 1 S 
MOMENTA  METHOD  PROGRAMS  FOR  SCATTERING  FROM  A  PERFECTLY 
CONDUCTING  (P*C. )  BOR.  DETAILS  OF  THE  MOMENT- METHOD  AND 
REQUIRED  (BUT  UNMODIFIED)  SUBROUTINES  ARE  DOCUMENTED  IN: 

1.  J.  R.  MAUTZ  AND  E.  F.  HARRINGTON,  'H- FIELD,  E-FIELD, 
AND  COMBINED  FIELD  SOLUTIONS  FOR  BODIES  OF  REVOLUTION* , 
TECHNICAL  REPORT  TH-77-2,  (FEBRUARY  1977),  AND 
•COMPUTER  PROGRAMS  FOR  H= FIELD,  E-FIELD,  AND  COMBINED 
FIELD  SOLUTIONS  FOR  BODIES  OF  REVOLUTION',  TR-77-3,  (MAY 
1977),  DEPARTMENT  OF  ELECTRICAL  AND  COMPUTER  ENGINEERING, 
SYRACUSE  UNIVERSITY,  SYRACUSE,  NY  13210 

THE  BOF  SURFACE  IS  DESCRIBED  BY  A  LONGITUDINAL  COORDINATE 
T  AND  AN  AZIMUTHAL  COORDINATE  PHI,  THE  FAR  FIELD  POINTS 
ARE  DESCRIBED  BY  THETA  AND  PHI  SPHERICAL  COORDINATES.  ALL 
SURFACE  VECTORS  AND  DYADICS  USE  THE  BODY  COORDINATE  SYSTEM 
WHILE  THE  VARIOUS  LINEAR  POLARIZATIONS  OF  THE  INCIDENT  AND 
SCATTERED  WAVES  REFER  TO  THE  SPHERICAL  COORDINATE  SYSTEM. 

ROUTINES  CALLED: 

MATB  -  CONSTRUCTS  PART  OF  THE  I MPBDANCE  OR  ADMITTANCE 
MATRICES  FOR  THE  IBC  EQUATIONS  FROM  THE  PERFECT 
CONDUCTOR  MATRICES, 

YZ  -  CALCULATES  THE  PERFECT  CONDUCTOR  IMPEDANCE  AND 
ADMITTANCE  MATRICES  (SEE  REFERENCE  1) 

PLANE  -  COMPUTES  MEASUREMENT  MATRICES  FOR  THE  TRANSMITTER 
AND  RECEIVER  ANGLES  (SEE  REFERENCE  1) 

DECOMP-  COMPUTES  THE  MATRIX  DECOMPOSITION  OF  THE  IMPEDANCE 
OR  ADMITTANCE  MATRICES  (SEE  REFERENCE  1) 

SOLVE  -  SOLVES  THE  LINBAR  SYSTEM  OF  ALGEBRAIC  EQUATIONS 
FOR  A  GIVEN  SOURCE  VECTOR  (SEE  REFERENCE  1) 

FILE  DEFINITIONS:  5  -  INPUT  FILE  (ALL  INPUTS  FREB  FORMAT) 

6  -  OUTPUT  FILE 
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INPUT  VARIABLES : 

NM  -  TOTAL  NUMBER  OF  FOURIER  MODES  (0, 1 , 2  ,.  .  .  ,  NM-1 ) 

NP  -  NUMBER  OF  PROFILE  BODY  POINTS  (ODD  NUMBER,  G.  E.  5) 

NPHI  -  NUMBER  OF  POINTS  IN  GAUSSIAN  QUADRATURE  INTEGRATION 
BK  -  WAVENUMBER  (UNITS  CONSISTENT  WITH  RH  AND  ZH) 

TT  *  TRANSMITTER  THETA  ANGLE  (DEGREES) 

(NOTE:  TRANSMITTER  PHI  ANGLE  =  0.0) 

P  -  BODY  PHI  ANGLE  (DEGREES)  IDENTIFYING  IBB  PLANE  IN 
WHICH  THE  SURFACE  CURRENT  DENSITIES  ARE  EVALUATED. 
TR  -  RECEIVER  THETA  ANGLE  (DEGREES) 

PR  -  BECEIVER  PHI  ANGLE  (DEGREES) 

ALP  -  INTBGRAL  EQUATION  WEIGHTING  COEFFICIENT 

RH  -  E  CYLINDRICAL  COORDINATE  OF  BODY  POINTS 

ZH  -  Z  CYLINDRICAL  COORDINATE  OF  BODY  POINTS 

X  -  ABSCISSAS  OF  GAUSSIAN  QUADRATURE  INTEGRATION 
A  -  WEIGHTS  OF  GAUSSIAN  QUADRATURE  INTEGRATIONS 
ZS  -  NORMALIZED  SURFACE  IMPEDANCE  (SCALAR,  CONSTANT) 

OTHBR  MAIN  VARIABLES: 

N  -  NUMBER  OF  EXPANSION  FUNCTIONS 

N2  -  MATRIX  DIMENSION 

ZST  -  T-T  COMPONENT,  NORMALIZED  DYADIC  SURFACE  IMPEDANCE 
ZSP  -  PHI-PHI  COMPONENT,  NORM.  DYADIC  SURFACE  IMPEDANCE 
Y  -  NORMALIZED  ADMITTANCE  MATRIX  FOR  THE  P.C.  NFIE 

Z  -  NORMALIZED  IMPEDANCE  MATRIX  FOR  P.C.  EFIE 

YY  -  NORMALIZED  ADMITTANCE  MATRIX  FOR  IBC  MFIE 

ZZ  -  NORMALIZED  IMPEDANCE  MATRIX  FOR  IBC  EFIE 

BT  -  MEASUREMENT  MATRIX  FOR  TRANSMITTER  THETA  ANGLE 

RR  -  MEASUREMENT  MATRIX  FOR  RECEIVER  THETA  ANGLE 

B  -  SOURCE  VECTOR  (FROM  INCIDENT  FIELDS) 

C  -  SOLUTION  VECTOR 

KT  -  DENOTES  TRANSMITTER  POLARIZATION  ( 1 =THETA,  2=PHI) 

KR  -  DENOTES  RECEIVER  POLARIZATION  (1 =THETA ,  2=PHI) 

NHEC  -  DENOTES  TYPE  EQUATION  (1=MFIE,  2=EFIE,  3=CFIE) 

TJ  -  T** COMPONENT  OF  NORMALIZED  ELECTRIC  SURFACE  CURRENT 
PJ  -  PHI-COMPONENT  OF  NORM.  ELECTFIC  SURFACE  CURRENT 
TM  -  T-COMPONENT  OF  NORMALIZED  MAGNETIC  SURFACE  CURRENT 
PM  -  PHI-COMPONENT  OF  NORM.  MAGNETIC  SURFACE  CURRENT 
E  -  NORMALIZED  RCS  ( SI G HA/LAMBDA* *2) 

OUTPUTS: 

1.  ALL  INPUT  PARAMETERS 

2.  SAMPLE  OUTPUT  OF  THE  EFIE,  MFIE  AND  CFIB  MATRICES, 

AND  SAMPLE  OUTPUT  FROM  PLANE,  DECOMP  AND  SOLVE. 

3.  FOR  EACH  INTEGRAL  EQUATION  AND  EACH  TRANSMITTER 
POLARIZATION,  THE  FOLLOWING  AFE  PRINTED: 

A.  TJ,  PJ,  TM,  PM  AT  N  DISCRETE  LOCATIONS  ALONG  THE 
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T  COORDINATE  (MEASURED  IN  THE  PHI=P  PLANE) 

B.  NORMALIZED  PCS  FOR  BOTH  RECEIVER  POLS. 

THE  MEANING  OF  'NORMALIZED'  IS  AS  FOLLOWS..  THE  IMPEDANCE 
(ADMITTANCE)  MATRICES  AND  THE  SORFACB  IMPEDANCE  ARE 
NORMALIZED  BY  DIVISION  (MULTI PLICATION)  WITH  THE  FREE 
SPACE  NAVE  IMPEDANCE.  THE  ELECTRIC  AND  MAGNETIC  SURFACE 
CURRENT  DENSITY  COMPONENTS  ARE  NORMALIZED  BY  DIVISION  WITH 
THE  MAVENUMBER.  IN  THE  PBOGRAM,  FH  AND  ZH  ARE  NORMALIZED 
BY  MULTIPLICATION  WITH  THE  WAVENUMBER.  THE  RCS  IS 
NORMALIZED  BY  DIVILION  OF  A  SQUARE  WAVELENGTH. 

THIS  VERSION  WAS  WRITTEN  BY  J.  ROGERS  OF  THE  RADAR  SYSTEMS 
GROUP,  MIT  LINCOLN  LABORATORY,  LEXINGTON,  MASS.  02173 
NOVEMBER  1983. 

DISCLAIMER:  THIS  PROGFAM  IS  EXPERIMENTAL  AND  WAS  WRITTEN  TO 

ILLUSTRATE  CERTAIN  PROPERTIES  OF  IBC  INTEGRAL  EQUATIONS. 

NO  EFFORT  WAS  MADE  TO  OPTIMIZE  COMPUTER  SPEED  OR  STORAGE. 
THE  BASIC  CODING  CAN  BE  VALIDATED  FOB  THE  P.C.  CASE  BY 
COMPARING  TO  THE  SAMPLE  PROBLEM  OF  REFERENCE  1.  FOR  THE 
IBC  CASE,  CHECKS  CAN  BE  MADE  WITH  THE  NUMERICAL  EXAMPLES 
OF  THE  TECHNICAL  REPORT  FOR  WHICH  THIS  LISTING  IS  PART  OF. 

COMPLEX  U,TJ  (120)  ,  P J  (120)  ,  E  (1 2)  ,SN,SNR,RT  (240)  ,RR  (240) 
COMPLEX  Y  (1600)  ,Z  (1600)  ,B(40)  ,C  (40)  , U1  ,CON JG 
COMPLEX  1H  (120)  , PH  (120)  ,  YY  (1600)  ,  ZZ  (1600) 

COMPLEX  ZS,ZST(20) , ZSP (20) 

DIMENSION  RH  (43)  ,ZH  (43)  ,X(20)  ,A  (20)  ,THT  (3)  ,THR  (3)  ,R2(20) 
DIMENSION  IPS  (40) 

READ (5,*) NM,NP, NPHI 

READ  (5,*) BK,TT, P,TB,PR, ALP 

READ  (5,*)  (RH  (I)  ,1=1  ,NP) 

READ  (5,*)  (ZH  (I)  ,1=1,  NP) 

READ  (5,*)  (X (K) ,K=1 , NPHI) 

READ  (5,*)  (A  (K)  ,  K=1 ,  NPHI) 

READ  (5,*)  ZS 

WRITE  (6,49)  Nil, NP, NPHI 
WRITE (6, 48)  BK,TT,P,TR,PR,ALP 
WRITE  (6,45)  (RH  (I)  ,1=1 , NP) 

WRITE  (6,44)  (ZH  (I)  ,1  =  1 ,  NP) 

WRITE (6,47)  (X  (K)  ,K=1,NPHI) 

WRITE  (6,43)  (A  (K)  , K=1 ,  NPHI) 

WRITE  (6,100)  ZS 

N2=NP-3 
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N=N2/2 

N3=N2+N 

N4=N2*N2 

N6=6*N 

0=  (0.  ,  1-  ) 

PI*3«  141592653 
P4=.0625/PI**3 
P8=PI/180. 

ET»«120.0*PI 

THT(1)=TT*P8 

P=P*P8 

THH  (1)  =TB*P8 
PB=PR*P8 
C 

DO  101  1=1, N 
ZSP(I)  =ZS 
ZST (I) =ZS 
101  CONTINUE 

DO  42  J*1 ,BP 
RH(J)=BK*BH  (J) 

ZH  (J)  =BK*ZH(J) 

42  CONTINUE 
DO  17  J*1,N 
U2  (J)  *1,/HH  (2*J+1) 

17  CONTINUE 

DO  54  J=1,N6 
TJ(J)=0- 
P  J  (J)  *0. 

TH(J)»0. 

PM (J)  =0. 

54  CONTINUE 

DO  55  J=1,12 
B  (J)  =0- 

55  CONTINUE 
KBITS  (6,9) 

C  LOOP  THROUGH  ILL  FOURIER  BODES 

DO  41  K*1 , NH 
NN*K-1 
PN=NN*P 
CS=COS (PN) 

SN*2,  •SIN  (PN)  *U 

PN=HN*PB 

CSB»COS(PN) 

SNB«2.*SIN(PH)*U 
IF  (BN,  EQ,  0)60  TO  56 
CS«2,  0*CS 
CSH»2.0*CSE 

56  LiNB-0 
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CALCULATE  YY  AND  ZZ  FOR  THIS  MODE 
CALL  YZ(NN,NP,NPHI,RH,ZH,X, A,Y,Z) 

00  130  1=1, N4 
YY(I)=Y{I) 

ZZ  (I)  =Z  (I) 

130  CONTINUE 
CALL  HATB  (N2,  ZSP,ZST,Z) 

CALL  MATB  (N2,  ZSP,ZST,Y) 

DO  131  1=1, N4 

YY  (I)  =  YY  (I)  +Z(I) 

ZZ  (I)  =ZZ  (I)  +Y  (I) 

131  CONTINUE 

C  LOOP  THROUGH  THE  THREE  IBC  EO'S 

DO  58  NHEC= 1,3 
GO  TO  (61,62,63)  ,NHEC 

61  DO  105  JJ=1 ,N4 

Y(JJ)=YY(JJ) 

105  CONTINUE 

WHITE  (6, 8 )  Y  (1)  ,  Y  (2) 

GO  TO  59 

62  00  103  JJ=1,N4 

Y  (JJ)=ZZ  (JJ) 

103  CONTINUE 

HEITE(6,8)Y  (1)  ,Y(2) 

GO  TO  59 

63  WHITE  (6,8)  YY  (1)  ,YY(2)  ,ZZ  (1)  ,ZZ  (2) 

DO  66  J=1,N4 

Y  (J)  =YY  (J)  +ALP*ZZ  (J) 

66  CONTINUE 

C  CALCULATE  HATRIZ  DECOMPOSITION 

59  CALL  DECCHP (M2, IPS, Y) 

WRITE  (6,8)  Y  (1)  ,  Y  (2) 

IF  (LAME.  NE.0)  GO  TO  57 
LANE=1 

C  CALCULATE  MEASUREMENT  MATRICES 

CALL  PLANE (NN,N,  1,THT,RT) 

CALL  PLANE (NN,N,1 ,THR,RR) 

WRITE  (6,8)  RT  (1)  ,RT(2),RR(1)  ,RR  (2) 

C  CALCULATE  SOURCE  VECTORS 

57  DO  27  KT=1, 2 
L= 2*  (NHEC-1) +KT 
GOTO  (31,32,33,34,35,36) ,L 

C  MFIE  -  THBTA  TRANSMITTER  POL. 

31  DO  21  J=1,N 
B  (J)  =-RT  (J+N3) 

B  (J+N)  =-RT(J+N2) 

21  CONTINUE 
GO  TO  53 
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32  DO  22  J*1,N 
JN*J+N 

B  (J)  *-BT  ( JN) 

B(JH)~ BT(J) 

22  CONTINOE 
GO  TO  53 

BPIE  -  THETA  POL. 

33  DO  23  J*1,N 
B(J)  *RT  (D) 

JN*J+N 

B  (JH)  *-BT  (JN) 

23  CONTINOE 
GO  TO  53 

EPIE  -  PHI  POL. 

34  DO  24  J*1,N 

B  (J)  —  BT  ( J+N2) 

B  (J+N)=ET  (J*N3) 

24  CONTINOE 
GO  TO  53 

CFIE  -  THETA  POL. 

35  DC  25  J*1,N 

B  (J)  =-HT  ( J+N3)  ♦  ALP*BT  (J) 

JN=J+N 

B  (JN)  *-BT  (J+N2) -ALP*BT  (JN) 

25  CONTINOE 
GO  TO  53 

CFIB  -  PHI  POL. 

36  DO  26  0*1 ,N 
J  N=  J  ♦  N 

B  ( J)  —BT  (  JN)  -AL P*BT  (  J+N2) 

B  (JN)  *-ET  (J)  ♦ALP*BT  (J+N3) 

26  CONTINOE 

SOLVE  HATBIX  BQOATION 

53  CALL  SOLVE (N2*  IPSVY#B,C) 

WHITE  (6#8)C(1)  ,C(2) 

CALCOLATE  NOB MALI ZED  ELECTRIC 
S DEFACE  COBBENT  DENSITY 

J1*(L-1)*N 
GG  TO  (11*12)  fKT 

11  DO  13  J*1 »N 
J2*J+J1 

TJ  (J2)  *TJ  (J2)  +C  (J)*CS 
PJ  (J2)  *PJ  (J 2)  +C  (J+N)  *SN 
13  CONTINOE 
GO  TO  14 

12  DO  15  J*1 , N 
J2*J^J1 
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T  J  (J2)  =TJ  (J2)  +C  (J)  *SN 
PJ  (J2)  =PJ  (J2)  +C  (J+N)  *CS 

15  CONTINUE 

C  CALCULATE  NOB  HA LI ZED  BCS 

14  DO  16  KR=1,2 
1=  (KT-1)  ♦2+KB 
KB=4*  (NHEC-1)  +L 
01=0. 

GO  TO  (71,72,73,74)  ,L 

THETA  -  THETA 

(TEANSHIT  POL.  -  RECEIVE  POL.  ) 

71  DO  75  J=1,N 
NPJ=N+ J 

01=01+ (BB  (J)  +  ZST(J)*RR  (N3+J))*C  (J) 

01=01+  (BB  (NPJ)-ZSP(J)*BR(N2+J))*C  (NPJ) 

75  CONTINOE 
E  (KE)=E(KE)  +01*CSB 
GO  TO  16 

C  THETA  «  PHI 

72  DO  76  J=1,N 
NPJ=N+ J 

01=01+  (BE  (J+N2)-ZST  (J)  *HR  (NPJ)  )  •C  (J) 

01=01+  (RF.  (N3+ J)  +ZSP  (J)  +RE  (J)  )  *C  (NPJ) 

76  CONTINOE 

B  (KB)  =B  (KB)  +Ol*SNR 
GO  TO  16 

C  PHI  -  THETA 

73  DO  77  J=1,N 
NPJ=N+J 

01=01+  (RE  (J)  +ZST  (J)  *RR  (N3+J)  )  *C  (J) 

01=01+  (RB  (NPJ)-ZSP(J)  *RR  (N2+J)  )*C  (NPJ) 

77  CONTINOE 

E(KE)  =E(KE)  +01*SNR 
GO  TO  16 

C  PHI  -  PHI 

74  DO  78  J=1,N 
NPJ=N+J 

01=01+  (RR  (J+N2J-ZST  (J)*RR  (NPJ))*C  (J) 

01=01+  (EB  (N3+J)  +ZSP  (J)  *RR  (J) )  *C  (NPJ) 

78  CONTINOE 

B  (KE)  =E(KE)  +  01*CSB 

16  CONTINOE 
27  CONTINOE 
58  CONTINOE 
41  CONTINOE 

C  WRITE  OOTPOTS 

DO  28  NHEC=1,3 
DO  29  KT=1, 2 
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WFITE  (6, 1 8)  NHEC, KT 
WRITE  (6,19) 

Jl=N*(2*NHEC*KT-3) 

DO  37  J=1  ,  N 
J2=J+J1 

T  J  (J2J  =TJ  (J2)  *R2(J) 

PJ  (J2)  =PJ  (J2)*E2(J) 

WpITE  (6,38)  TJ  (J2)  ,PJ  (J2) 

37  CONTINUE 
WRITE  (6, 119) 

DO  137  J=1#N 
J2*J+J1 

TM(J2)=-ETA*ZSP(J)*PJ  (J 2) 

PH  (J 2)  =ETA*ZST(J)*TJ(J2) 

WRITE (6, 38)  TM  (J2)  ,PH  (J2) 

137  CONTINUE 

DO  30  KR= 1,2 
Jl=4*NH£C+2*KT+KR-6 
SIG-P4*E(J1)*C0NJG(E  (J1)) 

WRITE (6,10) NHEC, KT,KR,SIG 
30  CONTINUE 
29  CONTINUE 
28  CONTINUE 
STOP 

FORMAT  ST ATEMENTS 

9  FGRMAT  ('0  SAMPLE  OUTPUT  FROM  SUBROUTINES*) 

8  FORMAT  (1X,4E14»  7) 

18  FORMAT  (»0NHEC=»  ,13,  «,  KT=»,13) 

19  FORMAT  ('  REAL  JT  IMAG  JT  REAL  JP  I  HAG  JP*) 

119  FORMAT ('  REAL  HT  IMAG  HT  REAL  HP  IMAG  HP') 

38  FORMAT  (1X,4E11<  4) 

10  FORMAT  (•  NHEC=',I3,‘,  KT=»,I3,',  KB=',I3, 

1»,  SIGMA/  (LAMBDA)** 2=',  Ell.  4) 

49  FORMAT  ('  NM  NP  NPHI  »/1X»2l3,I4) 

48  FGRMAT(7X,»BK',12I, *TT* , 13X,' P* ,1 2X, (TB • , 12X, • PV(/1X, 
15E14.  7/7X,  'ALP’/IX,  E14.  7) 

45  FORMAT  ( '  RH '/ (1  X,8F8.  4)  ) 

44  FORMAT  (*  ZH  •/ (1 X,  8F8.  4)  ) 

47  FORMAT (*  X  V(1X,5E14. 7)  ) 

43  FORMAT  (•  A'/(1X,5E14«  7)  ) 

100  FORMAT  (IX,*  ZS  =  ',2F10.4) 

END 

SUBROUTINE  MATB  -  THIS  ROUTINE  CALCULATES  TBE  MATRIX  (B) 
WHICH  IS  PART  OF  THE  NORMALIZED  IMPEDANCE  (ADMITTANCE) 
MATRIX  OF  THE  IBC  EFXE  (HFIE)  WHEN  THE  SURFACE  IMPEDANCE 
IS  NON- ZERO: 
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<B>  (I) 


1.0/BT1  (N)  (A)  (N) 


VHBBE  (A)  IS  1  P.C.  ADMITTANCE  (IMPEDANCE)  MATRIX,  (M) 
EEPBESENTS  THE  EXPANSION  COEFFICIENTS  OF  THE  MAGNETIC 
SOBFACE  CURRENT  DENSITY,  (I)  REPRESENTS  THOSE  OF  THE 
ELECTRIC  CURRENT,  (N)  BBPBESENTS  THE  CROSS  PRODUCT  OF  THE 
OUTNARD  UNIT  SURFACE  NORMAL  AND  ETA  IS  THB  FREE  SPACE  RAVE 
IMPEDANCE. 


SUBROUTINE  HATB  (N2, ZSP, ZST, X) 
COMPLEX  X(1),ZSP(1),  ZST  (1)  ,  CTEMP 
N=N2/2 
N2N=N2*N 


I  DENOTES  RON,  J  DENOTES  COLUMN. 
II,  12,  13  AND  14  ARE  IN  THE  TT, 
PT,  TP  AND  PP  SUBMATRICES  OF  X 
RESPECTIVELY. 

DO  1C  J=1 ,N 
11=  (J=1)*N2+1 
I2=I1*N 
I3*I1+N2N 
I4=I3+N 

DO  20  1=1, N 
CTEMP=X(I4) 

X  (14)  =X  (II)  *ZSP  (J) 

X(I1)=CTEHP*ZST  (J) 

CTEHP=I  (13) 

X(I3)=-X(I2)*ZSP(J) 

X  (12)  =-CTEMP*ZST  (J) 

11*11*1 

12=12*1 

13=13*1 

14=14*1 

20  CONTINUE 

10  CONTINUE 
RETURN 
END 
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GLOSSARY  OF  SYMBOLS 
(Order  of  Appearance) 


■v 

-% 


IBC 

MFIE 

EFIE 

CFIE 

BOR 

2s 

(fiSft1) 

S 

ej“t 

^vol 

(£8,fts) 

(J.ft) 

A 

n 

x 

(fi+ ,ft+) 


r 

n 

• 

Is 

y 

« 


Impedance  boundary  condition  (Leontovlch) 

IBC  magnetic  field  integral  equation 
IBC  electric  field  integral  equation 
IBC  combined  field  integral  equation 
Body  of  revolution 

Dyadic  normalized  surface  impedance  (double  arrow 
denotes  dyadic) 

Incident  electromagnetic  fields  (arrow  denotes  vector) 

Closed  surface  of  imperfectly  conducting  body 

Sinsusoidal  time  dependance  (phasor  notation) 

Volume  electric  current  density 

Scattered  electromagnetic  fields 

Equivalent  electric  (J)  and  magnetic  (ft)  surface 
current  densities 

Outward  unit  normal  to  surface  (symbol  over  n  denotes 
unit  vector) 

Vector  cross  product  in  equations 

Total  electromagnetic  fields  just  outside  (+)  or 
just  inside  (-)  equivalent  current  surface 

Field  point  vector  (source  point  when  primed) 

Free  space  wave  impedance 

Vector  scalar  product  in  equations 

Principal  value  surface  integration 

Permeability  of  medium  surrounding  body  (free  space) 

Free  space  Green’s  function 


r 


v 


e 

ir 

k 

I  I 

le.  lh 

a 

* 

t 

f 

ft 

[  1 
[V] 

[Z] 

[I] 

[U] 

[Y] 

[N] 

[M] 


GLOSSARY  OF  SYMBOLS  (Cont'd) 

Del  operator  (operates  on  source  point  coordinate  when 
primed)  used  for  gradient  (v*)»  divergence  (V*A)  or 
curl  (VxA) 

Permittivity  of  medium  surrounding  body 
Pi  -  3.14159... 

Wavenumber 

Absolute  value  (length) 

Generalized  integral  equation  weighting  coefficient 

Perfect  conductor  integro-dif f erential  operators 

Scalar,  constant  integral  equation  weighting 
coefficient 

Cylindrical  surface  azimuthal  coordinate 
Surface  coordinate  along  BOR  generating  curve 
Triangular  expansion  function 
Moment -method  testing  function 
Matrix 

Source  from  incident  electric  field  vector 

Normalized  impedance  matrix  for  perfect  conductor  EFIE 
(not  to  be  confused  with  surface  impedance) 

Electric  current  expansion  coefficients  (unknown) 

Source  from  magnetic  field  vector 

Normalized  admittance  matrix  for  perfect  conductor  MFIE 

A 

Matri/  representing  n  x  vector  operation 
Magnetic  current  expansion  coefficients 
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GLOSSARY  OF  SYMBOLS  (Cont'd) 


[ZZ] 

[YY] 

A 

V 

+• 

F 

[R] 

p 

[S] 
SEG 
Zs 
ka 

Re(  ) 
Im(  ) 
cond 
o 

a 

X 

0 

ME 


Generalized  impedance  matrix  for  the  IBC  EFIE 
Generalized  admittance  matrix  for  the  IBC  MFIE 
Electric  vector  potential 
Scalar  potential 

Magnetic  vector  potential 
Right-handed  spherical  coordinate  system 
Measurement  matrix 
Radial  cylindrical  coordinate 
Generalized  measurement  matrix 

Moment-method  segmentation  (segments /wavelength) 
Scalar,  constant  normalized  surface  impedance 
Sphere  wavenumber  -  radius  product 
Real  part  of  a  complex  quantity 
Imaginary  part  of  a  complex  quantity 
Approximate  matrix  condition  number 
Differential  scattering  cross  section 
Sphere  radius 
Wavelength 

Bistatic  scattering  angle 

Bistatic  cross  section  mean  error  (percent) 
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